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Abstract 



. The problem of positive Kolmogorov-Sinai entropy of the Chirikov-Standard map T\f : {x,y) i— > (2x — 

' y + \f{x),x) with f{x) = sin(a;) with respect to the invariant Lebesgue measure on the two-dimensional is 

(-H , open. In 1999, we believed to have a proof that the entropy can be bounded below by log(A/2) — C(A) with 

■ ■ C(A) = arcsinh(l/A) + log(2/V3) and that for A > Ao = (8/(6 - 3\/3))^''^ = 3.1547..., the entropy of TAsin 

should be positive. This approach was based on an idea of M. Herman using subharmonic estimates. 



1 The entropy problem of the Standard map 

The Chirikov- Taylor Standard map 

x\ f 2x~y + \f{x) 



(N 
> 

00 ; ^ \ V J \ X 

\ with J{x) — sin(a;) and real parameter A is a measure preserving difFeomorphisms on the two-dimensional torus 

■ = R^/(27rZ)^. It is probably the most famous example of a symplectic twist map. This map appeared in 
1960 in the context of electron dynamics in microtrons (see 1201) • It was first numerically studied by Taylor in 
1968 and Chirikov in 1969 (these independent studies are unpublished but see [371 E|) and also known under 
the name "kicked rotator" and describing ground states of the Frenkel-Kontorova model |79[[71. It is often used 
to illustrate or motivate various more general mathematical theorems in smooth dynamical systems (i.e. |131p 

^ \ the calculus of variations (i.e. [2111211) perturbation theory or renormalisation group techniques (for example to 

■ understand the break up of invariant tori [SHlfTSS, 136 ). While it is known that for A ^ 0, the map Txf is non- 
r> \ integrable, has positive topological entropy and horseshoes (i.e. [SEIEIEI), the question, whether hyperbohcity 

■ can hold on a set of positive Lebesgue measure stays open. Similarly, while many Lebesgue measure preserving 
diffeomorphisms on the torus are known to be non-ergodic with positive topological entropy (i.e. |140[ I157] l. 
it is not known whether positive metric entropy is dense in the C°° topology. The issue of the positivity of 
the Lyapunov exponents on some set of positive Lebesgue measure for Hamiltonian systems has been addressed 
at various places or reviews [H [133 [HHl [SB ISl [SI [EH [IHl 121111 1S2 GSSl S 11321 ■ According 
to |113l 15^ . the particular mathematical problem of positive entropy of the Chirikov Standard map had been 
promoted in the early 80'ies by Sinai. The textbook |131| states on p. 144 a conjecture (H2) that the entropy 
of the Chirikov Standard map is positive for all A > and that the entropy grows to infinity for A — > oo. 

Although numerical experiments show a very clear lower bound log(|A|/2) for the entropy of the Chirikov 
Standard map (e.g. |37lll9lfll5j ). it was not known even for one single value of A, whether the metric entropy 
with respect to the invariant Lebesgue measure can be positive. The second part of the (H2) conjecture of Sinai. 



Conjecture 1.1 



The Kolmogorov-Sinai entropy /i(T\ sin) of the Chirikov-Standard map T\s[ 
with respect to the invariant Lebesgue measure is hounded below by log(|A|/2) - 
C(A) where C{\) = arcsinh(l/A) + log(2/V3). 



'Department of Mathematics, University of Texas, Austin, TX 78712, USA (address of 1999) 
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Numerical experiments suggest that one could probably even hope to get rid of the C(A)-term in the case 
of the Standard map. Similar statements should hold for more general Standard maps and expect explicit 
lower bounds log(|mA|/2) — C(mA) if Asin(a;) is replaced by Ex + Asin(mx), where the integer i? e Z allows 
additionally to tune the homotopy type of the map. This property of positive metric entropy should be stable 
in the real analytic category and dense in the C'^(T) topology: 



Conjecture 1.2 



There is in C*'(T) a C'^ -dense set of real-analytic Standard maps / : i— s- R 
for which fJ^iTf) > 0. Let f be in this dense set. In every Banach space of 
realanalytic maps, in which f is, there is an open neighborhood of realanalytic 
maps g, for which fJ-{Tg) > also. 



The conjectured sensitive dependence on initial conditions could be persistent in the realanalytic category 
and should be obtained by realanalytic, C°(T)-small perturbations of integrable maps. The stability of pos- 
itive metric entropy with respect to realanalytic perturbations of the map would make the result physically 
relevant. Other results in Hamiltonian dynamics that have both this stability and which deal with orbits 
forming a set of positive probability is the theory of Anosov maps (see i.e. |81l 11521 IGSp (an example is 
{x,y) I— > {4:X + Xsin{x) — y,x) for small A) or KAM perturbation theory (see i.e. |l()5llK7ll(r7] 1 (which applies for 
example near the integrable SBKP map ( 139, 12 ) {x,y) i-^ (2a; — ?/-|-4arg(l-|-Aexp(— It is not known 
whether there are open sets of realanalytic Hamiltonian maps and flows for which there is quasiperiodic motion 
on a set of positive Lebesgue measure and simultaneously a conjugation to a Markov chain on a different set of 
positive Lebesgue measure. ^ 

One could ask similar questions to higher-dimensional symplectic maps {x,y) i— > {Ex ^ y + Xf{x),x) on T^'', 
where i? is a constant symmetric matrix in GL{d, Z) and / is a vector valued, real-analytic function on the 
torus. Examples are Froeschle maps, where Ex = 2x, fi{x) — A^ sin(a;i) -I- smi^^Xj) or classes of nearest 



neighbor coupled map lattices [HSI on T^*^, where Ex = 2x, fi{x) — Asin(a;i) -f e(sin(xi+i 
and XiJ^d = Xi. All averaged Lyapunov exponents should be nonzero for large enough A. 



' sin(a;i 



i)) 



As for any diffeomorphism T on a compact manifold M leaving invariant a smooth measure, the entropy 
is by the Pesin formula i92j equal to the sum of the positive integrated Lyapunov exponents. In the two- 
dimensional case, the entropy is lim„^oon^^ /^^ log | |(iT"(x, | dxdy. Since an estimate of the norm of the 
Jacobean of the Chirikov-Standard map T\ sin shows that the integrated Lyapunov exponent is bounded above 
by /j,2 log||dTAsin(a;,?/)|| dxdy < log(A/2) -f c(A) with c(A) = 0(1/A), Conjecture 11.11 provides for large A a 
rather accurate estimate for the actual values of the entropy of the Chirikov Standard map and shows that the 
measure of the Pesin region, the set on which the Lyapunov exponents is positive goes to 1, when |A| oo. 

By Pesin theory jll7l 1941 11181 , positive metric entropy (and in higher dimensions the non- vanishing of 
all Lyapunov exponents on a set of positive measure) would imply the existence of invariant sets of positive 
Lebesgue measure on which T\ is ergodic and on which some iterate of T\ is mixing and actually measure- 
theoretically conjugated to a Bernoulli shift. Other consequences are the density of periodic orbits in the Pesin 
region and shadowing properties (see [H^IH5] ). 



An obstacle for proving positive entropy estimates for the Standard map is Donskaya's observation that elliptic 
islands can exist for arbitrary large A ^1 1181 1861 021 ^^'^ which make it hard to find invariant cone bundles 
|148[I149] in the tangent bundle. In [32] it has been shown using renormalisation of coordinates near homoclinic 
tangencies 11041 that for a residual set of large parameters A one has Standard maps TAsin which are nonergodic. 



^There is a possibility to achieve such mixed behavior on the union of two closed symplectic manifolds, where the map is Anosov 
on one component and integrable on the other. For smooth, not realanalytic examples in the literature that are unstable under 
perturbations of the map. 
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Chirikov had already expressed concerns that the 
measure of the stable elliptic component might be 
close to 1 p. 333). It has been often asked 
whether an area-preserving monotone twist map 
has a dense set of elliptic islands in general (see 
e.g. |52l It is known that a Baire generic 

symplectic non-Anosov C^-diffcomorphism has a 
dense set of elliptic periodic orbits |1U9| . Some re- 
search towards avoiding elliptic islands by smooth 
ergodic perturbations of the Chirikov Standard 
map has been done in |132| . It is possible to get 
positive entropy by a smooth perturbation of 
the map [55]. Whether this is possible with C°° 
perturbations is not known |55| . 



The question whether a dense set in Diff^(Af) of measure preserving diffeomorphisms on a manifold AI has 
positive metric entropy has been asked in p)9; . 

It has been conjectured that there exists a set of parameters A with full density at oo for which the Chirikov 
Standard map has no elliptic islands ^Hl (see also |32p. 

Our proof-attempts of positive metric entropy depended on spectral and complex analytic techniques as well as 
the determinant theory of finite von Neumann algebras avoiding the ergodicity question. We hoped to obtain 
many non-ergodic realanalytic Standard maps with positive metric entropy. The picture above for example 
shows plots of some orbits of the Chirikov Standard map in the case A = 3.4, where a linearly stable fixed point 
(1/2, 1/2) coexists with a region with positive metric entropy. 

There are Standard maps T^sm arbitrarily close to this map 13.48111 for which the Birkhoff normal form at 
(1/2, 1/2) is such that the fixed point is surrounded by invariant KAM curves of positive Lebesgue measure and 
for which we still have ergodic components of positive Lebesgue measure. Together with we know that 
there exist for large A a open dense set of parameters for the Chirikov Standard map which lead to nonergodicity 
and positive Kolmogorov-Sinai entropy. 

Consider the matrix-valued map 

Together with a Lebesgue measure-preserving dynamical system T : ^ T^, it defines a cocycle {x,y,n) 1-^ 
A/ t(^' y) ~ ^E,xfiT"'~^{x, y)) o • • • o AE,\f{x, y) for which the integrated Lyapunov exponent 

K^E.\f.T) = lim - / logP^ w -r(x,?/)|| 
is defined. For general information on Lyapunov exponents, see for example |125l I84| . 

Pesin's formula links the entropy of Tf with the Lyapunov exponent of dTf and so with the Lyapunov exponent 
of the cocycle because AE=2,df,Tf agrees with the Jacobean dTf of the map Tj. 



2 A nonselfadjoint spectral problem 

Consider the random operator 

{Lw{x,y))u„ = u„+i + Un-i + Xiw^^ exp(ia;„) -I- w exp(-ix„))u„ , (1) 
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on P{Z,C), where {xmVn) — T^^{x,y) and w e C is a complex parameter. For fixed {x,y) e T^, this is 
a bounded operator on Z^(Z,C). The term 'random' is used because {x,y) L{x,y) is an operator valued 
random variable. Such operators are in a von Neumann algebra with finite trace tr{K) — Jip2 [^(a;, y)]oo dxdy 
(see e.g. [THZIEHIEI])- 

Similar as in the case, when the complex parameter is the energy, the Lyapunov exponent ^{w) of the transfer co- 
cycle Ae{w) is a subharmonic function in w. As we will see below, the Thouless formula ^{w) — log(det(L^)) = 
tr(log \Lw\) is still true for the in general nonselfadjoint operator L^. 
The problem is to compute det(L^). The entropy of the Standard map is logdetLu,=i- 

The average value of the Lyapunov exponent ^.{w) on the circle |ti;| — 1 can be estimated with the subharmonic- 
ity argument of Herman 56 or with a Jensen formula as used by Sorets-Spencer in |133| . In higher dimensions, 
an adaption ,48, of ^133, to higher-dimensional cocycles would apply. 

A theorem of Lax "F? on continuous measure-preserving transformations allows to approximate the group X of 
measure preserving homeoniorphisms by finite groups yu of measure preserving transformations. If one looks at 
the value of the Lyapunov exponent of the matrix-valued map AE.\sin,T on J^fe, (where evaluating the Lyapunov 
exponent is a reliable finite dimensional integration), one can observe numerically that for every transformation 
T, one can find a transformation T such that ^{AE,\sin,T) + l^-i-^E Asin f) — 21og(A/2). This relation can be 
understood in terms of determinants. 

By looking for a proof of this relation, we were motivated by Aubry duality which can be defined if T is in 
such a finite group 3^^. Aubry duality is an involutive transformation L L, which preserves the density 
of states such that log(det(iA)) ~ log(det((A/2)L(4/;)j))) = log(A/2) + logdet(L(4/A)). In the Mathieu case 
T{x,y) = {x + y,y), where Lx — L\, this leads with log(det(L4/;^)) > to the estimate logdet(LA) > log(A/2) 
which is expected to be true also if T is the Standard map. Here, we will show that there exists for each 
transformation T and corresponding operator L = Lt,\ a different transformation S and operator L — Ls^\ 
such that log(det(LL)) > 21og(A/2). 

How is the operator L obtained? The group Ra ■ {x, y) ^ (x + a, y) of translations on the torus acts by conjuga- 
tion on the whole group y of measure preserving transformations. The orbits of this action are T^q, — Ra^R—a- 
The operator L will be obtained by changing T to one of the conjugates Ta ■ The value of the Lyapunov exponent 
w = exp(iQ!) t—^ fi{AE,\f.Ta) extends to a subharmonic function from the circle \w\ = 1 to the entire complex 
plane C. If one considers the complex parameter w as a, 'spectral parameter', the Riesz measure dk used to 
represent the Lyapunov exponent plays the role of the density of states. It is a measure in the complex plane 
which has in our case in general positive area in the complex plane. The measure dk has its support in an 
annular neighborhood of the unit circle (and an atom at 0). 

We use the artificial spectral parameter w because with respect to z = re'^, there is no analyticity of 
[z" A" (z , y)]ij . The motivation is trying to generalize results for maps which extend analytically to polydiscs 
like for T{x, y) ~ (x + q, y), where both Aubry duality |S] as well as Herman's subharmonicity |56| work directly. 
While the map (te*^,e*^) i-^ (j^^i{2x-y+\sin{x)) ^ ^ix^ similar complexification attempt) is real-analytic in 

each of the variables t, x, y, it is definitely not analytic with respect to z = te" and functions which needed 
to be subharmonic are not. Our first attack on the Standard map using plurisubharmonicity was done in the 
spring 1988, now more than eleven years ago. The approach was the observation that the Standard map family 
written in Hamiltonian form {x,y) ^ {x -\- y + Asin(x),?/ -|- Asin(a;)) is induced on invariant tori of the single 
analytic map U : (z, w, u, v) i-^ (zwe^"", we^"", wwe"~^, we"~^) on C^, and where nonanalyticity manifests itself 
that the t/-invariant two-dimensional tori S\ = {{z,w,u,v) \ \z\ — \u\ — A/2, \w\ = \v\ — 1 z — u,w = v} 
on which U induces T\, are not distinguished boundaries of polydiscs. Later attempts were to Vlasov-Toda 
deform the cocycle with the aim to estimate the Lyapunov exponent on tori of isospectrally deformed operators 
(on which the Lyapunov exponent is constant) |71[I7(J| or to fix the dynamical system and to vary the cocycle |69|. 

Because w i— > gn{w,x,y) = w"'[A^{w, x,y)]ii is analytic, the r — \w\ dependence of the local Lyapunov 
exponent can be related with the average angular dependence of the argument: in polar coordinates w = re"^. 
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the Cauchy-Riemann differential equations 



4l- f arg(5„(w,x,y)) dxdy = -7- / log\gn{w,x,y)\ dxdy 
d(f> r J'j-2 dr 

are valid in any connected region of the w- resolvent set. In the complement, in the support of the w-spectrum, 
this formula can not be used and an integrated version, a Jensen formula replaces it. Besides the contribution of 
the radial change of the argument, there is an additional nonnegative subharmonic contribution to the Lyapunov 
exponent. 



3 More heuristics 

T , A ( \ ( E+^{z + z-^) -1 \ 
l^et AE,\cos[z,y) = \ Q I. 

In situations like the Mathieu case T{x,y) = {x + y,y), where T extends to an analytic map on D x T with 
respect to z = rexp(zx), Herman's subharmonicity 56 gives /i(^T) — log(A/2) > 0. For a general map T G y 
which does no more commute with i?^ : {x,y) 1— *■ (x + a,y) on T^, the function w = exp(iQ;) fJ.{AE,\sin,Tc,) 
with Ta = RaTR-a is no more constant. In average, we have /-j. ^(^E.Asin.Tc) da > log(A/2). Furthermore 
/^(^B.Asin.Ta) is bounded above by log(A/2) + c(A) with c(A) = 0(1/A). This means that for large A, we have 
fJ-{^E,Xsin,T^) > log(A/2) for a large set of a's. We used such a fact in a similar way in |^ . 

However, we don't know the specific value of the upper-continuous function a '—^ fi{AE,XBin,Ta) the point 
a = 0, which we are interested in. It could be zero a priori. Numerical computations however show that the 
minimum of a 1— > fJ-{AE^Xsin,Ta) is not farther away from the mean value than the maximum. There is a heuristic 
explanation which uses the Jensen formula in a sector and which will assume T £ X. This Jensen formula (see 
Section l^l is essentially an integrated Cauchy-Riemann differential equation and sharpens the subharmonicity 
tool. We also use the Lax approximation theorem in Section [TUl to explain the now following heuristics. 

The Lyapunov exponent of the matrix- valued map {x,y) G t—f _B(e*^) with 

\ A f \ ( E+^{z + z-^) -1 \ 
B{z) ^ zAe^XcosKz) ^ z\ / 

is the same as the Lyapunov exponent of the cocycle {x^y) i— > Ae,Xcos{x) because \z\ = 1. 
Define for (x, y) G and n e N the complex function 

w i-> 9n,T{w~^ exp{ix), y) = [Be^xcosAW^ exp(ia;), ?;)]ii 
which is analytic in the complex plane. 

Denote by Argj^^^ z2]i9n,T) the argument change of the analytic function gn.T on the line from zi to Z2. Unlike 
arg(g„,T), the argument change Argj^^ z2\^9n,T) is a uniquely defined number. 

Let Tk be a cyclic Lax cube exchange transformation approximating a homeoniorphism T £ X . The map Tfc is 
extends from to D x T by setting rfc(r exp(ia::), j/) — (r exp(ia::i), yi) with Tk{x,y) = 

In each sector, we could also take the z variable instead of the w variable because z ^ A"{z^y) is analytic in 
each sector As,tj = {w e C \ \w\ e [s,t],a.Tg{w) £ [27rj/fc, 27r(j + l)/fc) }. 
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Let r+(Q) be path connecting se'" and ie*" on 
a radial line segment, T~ the reversed path. Let 
7+(q!,/3) — {w,\w\ — s,arg(w) g be a 

circular arc contained in the circle 7s — {\w\ — 
s} and 77 (a,/?) the reversed path. The differ- 
ence between the argument change of the func- 
tion gn.xiw) = [■B£;,Asin,T("')]ii along the cir- 
cles 7 J and enclosing the annulus As,t — 
{|w| G [s,t]} is the sum of the k angular 
argument changes ^T:g^+i^2-.j/k,2^u+i)/k)'^9^-T) " 
Arg^^-(2xj7fc,27r(j+i)/fc)(5",T) plus the sum of the 
k radial argument changes Argp+(27rj7fc) (5n,T) — 
Argr-(27rO+i)/fc)(5n,T) with j ^l,...,k. 



The total sum of the radial plus angular argument changes is nonnegative because it is the sum of the indices 
of the roots of z t—^ gn,T inside the annulus As,t^ all of which are nonnegative. 

The sum of the k radial argument changes is bounded above by a 0(l/A)-term, because a positive value 
contributes positively to the Lyapunov exponent estimate and the Lyapunov exponent is bounded above by 
log(A/2) + c(A) with c(A) — 0(1/A). There is no reason, why the sum of the k radial argument changes should 
differ much from its negative. Numerical experiments strongly confirmed that they don't. If this sum of the 
radial argument changes 

k 

^Ar gr+(27rj7fc)(5n,T) " Argf- (2^(j + i)/fc) (ffn^r) (2) 

(which is averaged over and bounded uniformly in k) would vanish, one would have the same estimates as 
in the subharmonic case and 

""W ^og\9n,Tix,y,w)\ dxdy > n~'^ log |.g„_T(a;, y, 0)| dxdy = log(A/2) 
for all \w\ — 1. This would imply IjI.{Ae,\cos,t) > log(^/2). 

While the sum of the radial argument changes (O does not vanish for a general T E y, numerical data makes 
^ appear to be small. Indeed, we will see that for every transformation T E y, there is an other transformation 
Sa — RaTR-a such that the Lyapunov exponent of At plus the Lyapunov exponent of As beats 21og(A/2). 
This will imply that the fluctuations of 

^i,,{a) = n-^ [ logW {e'^"+^\y)]n\ dxdy 
around its mean /^/in(a) da > log(A/2) can not get too large. 

We have not been able to prove such a statement using complex analytic methods only. 

4 The spectrum of periodic difference operators 

This section revies some facts about not necessarily selfadjoint periodic Jacobi matrices and higher order dif- 
ference operators. Unlike in the selfadjoint case, when the spectrum is confined to the real axes, the topology 
of the spectrum can then be more interesting. All we need to know for our purposes is that the spectrum is 
contained in a nowhere dense set. 

While nonselfadjoint differential operators have been studied quite a bit by Russian mathematicians, nonselfad- 
joint higher order difference operators have appeared less frequently in the literature. But they do occur: these 
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operators appear naturally when studying the stability of stationary solutions in coupled nonharmonic oscilla- 
tors |144[l62j . For nonself adjoint Schrodinger differential operators, see 0^1 and references. 



Jacobi matrices. Three vectors a, 6, c € define periodic sequences a„ = a„+p,6„ = 6„+p,c„ = c„+p € 
Z°°(Z, C) which can be used to build a p-periodic Jacobi matrix 

{Lu)n = a„Mn+l + C„_iU„_i + 6„U„ 

which is a bounded operator on on P{Z, C). Denote by Lper the same operator L but acting on a different Hilbert 
space, namely the finite dimensional space of p-periodic sequences equipped with the norm (X]r=i l^ip)^^^. For 
w € C, (a parameter which has nothing to do with the complex parameter w in the operator ^ ), this operator 
is given by the matrix 

/ 5i ai • w~^Cp\ 
Ci 62 02 ■ ■ 



Lper(w) 



C2 • • 

• • Cp_2 O-p-l 

Vuia^ • c„_i bp / 



-i,p \j \j ^p- 

and satisfies Lp^r — Lper{l)- Denote by A(z, w) — det{Lper{w) ~ z) the characteristic polynomial of Lper{w) and 
by A(z) — A(z, 1) the characteristic polynomial of Lper- Comparing coefficients in w and using A(2;, 1) = A(z), 
one gets 

A(2, w) = A(z) — aw — cw~^ — b 
with a = 11^=1 % ~ 11^=1 '^i ^ ^ ~ Denote by r the shift ii„ 1-^ u„+i on ?^(Z, C). 

By the theorem of Cayley-Hamilton, we have A{Lper) = 0. We deduce from this that the operator A(L,tP) on 
Z^(Z, C) defined by the functional calculus (plug-in the operator L for z and the operator for w in A(z,w)) 
satisfies 

= A(iper) = A(i, ) = A(L) - ar?' - cr^f - b . 

This implies that the operator K :— A(i) = ar^ + c{tP)* -f- & on P{Z,C) is a Laplacian with space indepen- 
dent entries. Conjugated by a Fourier transform U : P(Z,C) — > L^(T, C), it is unitarily equivalent to the 
multiplication operator 

{UKU*)u{0) = (a e'P^ + b + c e-'P'^)u{9) X{9)u{0) 

on L^(Z,C) which has as the spectrum on an ellipse £ possibly being degenerated to an interval. Since 
A(L) = K, we know by the spectral theorem that a{A{L)) = A{(t{K)). The spectrum of L consists of p closed 
real curves. They are the zeros of z i-^- A{z) — when w is runs through the ellipse £. The density of states 
of L is the push-forward measure (A^^ o X)*dd of the Lebesgue measure dO on the circle T = R/(27rZ) under 
the multi-valued map — A^^ o A : T C. (See for a derivation in the case an — Cn — !)■ 

Remarks. 

1) Since the spectrum of Jacobi matrices is simple, the curves Zi{6) do not intersect transversely. However, by 
deforming along a one-parameter family of Jacobi matrices, curves could merge or separate. The discriminant 
{L I (Az,)'(0) = for some z € cr(i)} is the set of Jacobi matrices, where the topology of the spectrum can 
change under small parameter variations by merging or separation of spectral arcs. 

2) The Bloch variety of L is the set BViL) = {(z, w) G C* x C* | A(z) ^ aw + b + cw'^}, where C* is the 
Riemann sphere C U {00}. It is the set of all {z,w) such that there exists a nontrivial solution of Lu — zu 
satisfying (r^ujn — Wn+p — wu„. The spectrum of L is the intersection of BV{L) with {\w\ = 1}. 

3) A theorem of Polya (see 0) implies that the projection of the spectrum of a periodic (not necessarily self- 
adjoint) Jacobi matrix onto any real line in the complex plane has Lebesgue measure < 4. 

4) In the case when aj, bj, Cj are real and a„ — Cn and a = Y[j ctj = Ij the Jacobi matrix L is real and selfadjoint 
and K is the free Laplacian t + t* — 2 with spectrum [—2, 2] so that cr(i) = {z \ A{z) S [—2, 2]}. The spectrum 
is then a union of bands on the real line which can intersect only at the boundaries. 
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5) In the case when a — Y\- aj ~ 0,c — Y[j bj — 0, the spectrum is A~^(0), the set of roots of the characteristic 
polynomial of Lper confirming that in this case, L and Lper have the same spectrum. {L is then a countable 
direct sum of operators Lper-) 

Higher order diflference operators. For a p-periodic operator L — X]i=-r '^''''N the spectrum of Lp^r is more 
complex. The spectral bands curves can cross and for selfadjoint matrices, the bands can overlap. 

We consider the case p = fciV, where the operator can be rewritten as a Jacobi operator L = ar + ct* + on 
P(Z, C^), where Cj, 6j, Cj are N x N matrices satisfying aj^k — cij, bj+k — bj, Cjj^k = cj. 

Lemma 4.1 (Spectrum of periodic matrix-valued Jacobi matrices) The spectrum of the Jacobi opera- 
tor L with N X N -matrix valued coefficients is contained in the union of at most p — kN closed real curves. 

Proof. Unlike in the case = 1, the Bloch variety 

{(z, w) = A(z, w) = det(Lper(w) - = 0} C 

can in general no more given explicitly in the form f{w) = g{z). The spectrum is still the intersection of this 
complex variety with {\w\ = 1}, the parameter \w\ = 1 labeling Bloch waves. 

The intersection of the Bloch variety with = 1 is a set in the complex plane which is contained in a union 
of at most p — kN closed curves. 

Let Det denote the determinant function on the algebra of complex N x N matrices. We write 11^=1 — 
afc . . . ai, also in the noncommutative case. With a — Det 11^=1 o,j,c = Det 11^=1 '^j' ^e have 

A(z, w) = A{z) + aw^ + cw"^ + f{z, w) , 
where the mixed term f{z, w) is a polynomial in w, , z which is of degree < in it;, and of degree < p 

Examples. 

1) We are especially interested in higher order difference operators which arethe product L — L'^^^L'^^^ of two 
Jacobi matrices L'*^ = + t* + V'^'^K 

2) If a„, 6„, c„ are diagonal, the operator L is the direct sum of N one-dimensional Jacobi matrices. There are 
therefore operators with p — kN spectral curves. 

5 Thouless formula for difference operators 

This section reviews the Thouless formula for operators on the strip, that is Jacobi matrices with matrix-valued 
entries. We recall the proof in [HO] because the class of operators differs a tiny bit from the operators consid- 
ered in |Hn|: the off diagonal matrices are the identity to which some nilpotent part is added. An other change 
is that we do not assume selfadjointness of the operator. The proof of the Thouless formula however is the same. 

The Thouless formula ij.{Ae) — J log \E — E'\ dk{E') relates the density of states dk of a random Jacobi operator 
T + T* + 6 with the Lyapunov exponent ^{Ae) of the transfer cocycle (see |2S1E1)- This generalizes to operators 
r + T* -1- 6 on P(Z, C^), where tm„ = m„+i and b{x) is a selfadjoint N x N matrix, if ^(^b) is the average of 
the N largest Lyapunov exponent of the 2 A'' x 2Ai"-cocycle |Sn]. The Thouless formula extends also to the case 
of operators ar -I- (ar)* -f b, where a{x, y) are unitary |74| . 

The next lemma provides a generalization to operators of the type L — i'^^L^^) with Jacobi matrices 

(L('=)(a;('=\2/('=))u)„ = u„+i + + V;W(:rW, j/C^))^™ , 
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where the later are not necessarily selfadjoint operators on P{Z, C), 
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(3) 



c„(x) = V;i^^(x(i),?;(i)) + V;^^\(x(2),y(2)). The potentials vJi^\x'-J\y'-J^) are obtained from Tj e and func- 
tions y(^) : R by = V'^^^TJ'{x(j\y'^^^)). 



Each difference operator operator L(x) on P(Z, C) can be written as a Jacobi operator 

L = af + (cf )* + b, {Lu)n = flftUfi+i + Cft_2Ufi-i + 
on /^(Z,C^) with rUn — Ufi+i,fc(x) = c(rj^(a;(-^\ y^^^), r2^(a;(^\ y(^)))f and with matrix-valued entries 
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where n = 2ri. The transfer cocycle 
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if (L(x)u)„ — Eun for u G /^(Z, C^). The density of states of L defined as the functional / [/(_L(x)]oo dx. 

on C(C) and the density of states of L defined as the functional / ^ J!p2 Trace([/(Z(x)]oo) x both exist and they 
are the same probability measure in C, if Trace is the normalized trace on 2 x 2 matrices satisfying Trace(l) = 1. 

Remarks. 

1) For L = ^j^g^g 

are Jacobi operators, there are other transfer matrices like the product of the 
individual 2x2 transfer matrices of the operators X'*-'. These are different matrix- valued functions in E and 
the Thouless formula would look different in this case. 

2) Note that dk is normalized to be a probability measure. In some of the literature like for example in |8(Jj . 
the usual Trace for N x N matrices is taken which implies that the total mass of the density of states dk for 
operators L on the strip becomes N. 

Lemma 5.1 (Thouless formula for operators on the strip) Let dk be the density of states of the opera- 
tor Then 

f log \E ~ E'\ dk{E') = h{A^Ae)/2 , 
Jc 

where ^{a'^Ae)/'2 is the arithmetic mean of the two largest Lyapunov exponents of the 4x4 transfer cocycle 
Ae of the operator 

Proof. Small changes are needed to the known proofs |28l 1171 IHH] for the same statements for selfadjoint 
operators L, where a„ = c„ = l,rt G Z. Write x = (a:'^^^ x'-^-' , y^^^) for a point on T'*. 

(i) It is enough to prove the statement for periodic, (nonergodic) transformations T x S. 

Proof. Every T x 5* G 3^ x 3^ can be approximated by periodic transformations T„ x 5„ using Rohlin's lemma 
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|53 | (which implies that there exist measurable sets Yn whose Lebesgue measure on goes to 1 for n G oo so 
that (T X S')'=(x) = (T„ x Sn)^^) for fc = 0, . . . n - 1). 

The Thouless formula is an identity for subharmonic functions fxxs — gxxSi parameterized by T x S" G 3^ x 3^. 
Let dk{f) denote the Riesz measure of a subharmonic function / it satisfies dk{f) — A/ in the sense of 
distributions, where A is the Laplacian. It follows from the Avron-Simon lemma and the dominated convergence 
theorem that if r„ x S'n — > T x S" in the uniform topology for measure preserving transformations, then 
dk{fT„xs„) dk{fTxs) weakly. Also dk{gT„xS„) dk{gTxs) weakly holds because dk{gT) is a weak limit of 
measures dk{g!^'') — J2j where Ej''\x.) are the roots of i-^ [A^^(x)]ii and T i-^ dk{g^p^) 

is continuous from 3^ with the uniform topology to the space of measures on C with the weak topology. 
By assumption, we have rffc(/T„xS„) = dk{gT„xsJ- The refore dkifxxs) = dk{gTxs) which is A/txS = ^Qtxs- 
From this, JtxS ~ gxxs follows by Weyl's lemma jl21) . 

(ii) Next, it is sufficient to assume that E is outside the spectrum of a periodic operator L (which is now 
contained in a finite union of one-dimensional curves (Lemma (|4.1() '1'1 as well as that E is outside the spectrum 
of L'^^\ where L'^^) is the Jacobi matrix on Z(Z+,C) with zero Dirichlet boundary condition at n = 0. The 
spectrum of L'^^^ is a finite set of points. 

Proof. This is the Craig-Simon subharmonicity argument Taking balls Br{E) around a point E in the 

spectrum, the result holds for the over Br{E) averaged Lyapunov exponent. Passing to the limit r — > is 
possible because of subharmonicity. 

(iii) li\j{E) = YYk^iiE- \jk)-,i = 1, . . . , is an eigenvalue of the A^ x A^ matrix (where the 27V x 27V 
matrix Ae is written as a 2 x 2 matrix oi N x N matrices), then each root Xjk oi E i-^ gives rise to an 
eigenvector m^^'^' of a truncation L*^"' of L, where ak, bk, Ck are put to zero for fc < as well as for k > n. The 
measures dkn = (2n)~^ f. S{Xjk) converge for n ^ oo weakly to the density of states dk which is supported 
on finitely many curves. 

This is the Thouless- Avron-Simon lemma |23) which holds in general for random finite-difference operators 
which do not need to be selfadjoint. 

(iv) The fact that 

/ \og\E-E'\ dk4x,E') ^J2^og\E - A,,(x)| = log |Det([A'^]n(x))| 

converges to /z(A^A£;(x))/2 for almost all x e T'* in the limit n — !■ oo is seen in the same way as in the 
Appendix of IHOl: If m±(x) = a{x)u^{{T x S')x)u± (x) ^ are the matrix-valued Titchmarsh-Weyl functions, 
where I/tt^(x) = Eu^{yi) and {M^(x)}„g±z+ e P{±Z+ , M{N,C)), then 

(t1/2) / log|Det(m=^(x))| dx 

is the arithmetic mean /Lt(A^A£;(x)/2 of the two largest Lyapunov exponents of Ae- Furthermore, the maximal 
Lyapunov exponent averaged over is 

— / log ||to'''(x)|| dx = / log ||m~(x)|| dx . 

Because E is not an eigenvalue of L{x), there are two A^ x A^-matrices U, V such that the functions m^, m^"-* = 
au(°)(T X 5)(m^°^)""'^ satisfy m^^'' — m^U + ni^V . Here it'^°)(x) is an other solution of L{x)u = Eu satisfying a 
Dirichlet boundary condition at 0. The matrix U is invertible because else there would exists a 7^ with Ua = 
and m'^'^^a = rn^Va would be in P(Z+,C^) and contradict the fact that E was also chosen away from the 
spectrum of i*^"^. (In the selfadjoint case, this would already have been taken care of by choosing Im(£^) 7^ 0.) 
Writing 

ml") = {l + mtVU-\m-)-^)m-U 
and using mn\^) — [A^]ii(x), one concludes 

logdet(L) = / \og\E- E'\dk{E') 
Jc 
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lim / / log|£;-£;'| dfc„(x,^') 
lim / logDet([A|;(x)]n) dx 
i lim / log|Det(m(°)(x)„)| dx 

log|Det(TO"(x))| dx = ^i{A^Ae)/2 . 



fp4 



6 An upper bound for the Lyapunov exponent 

In this section we determine the constant C(A) in Conjectm'e and formulate a Proposition which outlines, 
how one can improve the constant with more effort. The reason to stick to the constant C(A) is that it allows 
explicit expressions. We also show in this section that C(A) can be replaced by c(A) = 0(A~^). 

Because ||A"|| < ||A||" for any matrix A, the number 

m{E,\) := / \og\\AE,xcos{x,y)\\dxdy 



is an upper bound for the Lyapunov exponent of the Standard map, rsp. the Lyapunov exponent of the cocycle 
Ae{w) for \w\ = 1. Define C(E, A) = m{E, A) - log(A/2) and C(A) = C(0, A). 



Lemma 6.1 C{E,\) = log(2/V3) + 0(1/A) for large A and C(0,A) > /or A > Ao = 2j2/(6-3V3) = 
3.1547.... 

Proof. With the norm |||A||| = max^^i 2 |^e,;|, where are the basis vectors, there is an explicit expression for 
an upper bound of the integral m{E, A). Following I6UI . define 



M{E, \) = \E + i + v/(S + i)2-A2|/V3 
where the square- root takes the solution with positive imaginary part. 

For matrices A = \ \ _ one has |||yl||| IVl + c2| > (\/3/2)||A||. Furthermore, with the analytic 



^-10^ 

function f(z) — z^iX + {E — A(z + z^-'^)/2)^) ~ X/2{z — a_)(z — a+) where |a_| < 1, one has using the Jensen 
formula 

\o%\\\AE,XcoAx)\\\dx = I- [ \og\l + {E-\coaix)f\dx 



fp 2 ^ fp 

1 
2 



log|/(z)| dx 

{\z\ = \ oxp(ia:)| = l} 



i(log(/(0))-log(|a_|)) 
log(^M(i?,A)). 



Therefore C{E,X) = \og{M{E, A)) - log(A/2) = log(2/ ^3) + log(| {E + i)/X+ y'iE + i)^/X^ - 1\) and C(0, A) = 



C{A) = log(l/A+^l/A2 + l) + \og{2/V3) = arcsinh(l/A) + log(2/V3). For A = 3.1547.. = 2y^2/(6 - 3\/3), 
one has C(A) = 0. 



The function A i— > C(A), having a positive derivative, is strictly monotone. 
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Lemma 6.2 We have 



[ \og\\dTE,xcos{x,y)\\dxdy<c{E,X):^ [ log + {E + Xcos{x)) dx = OiX'^) . 

Proof. We use Hadamard's determinant theorem, which tells that for a finite matrix L, one has logdet(i) < 
llcfjlb, where aj are the rows of the matrix L = (ai, • • • ,a„). (The relation follows from the fact that 
if A = (oi, . . . ,a„), then det(A)/ J^J^^ IWjlU = det(i3), where B = (ai/||ai||, . . . ,a„/||a„||). The determinant 
theorem is equivalent to det(i3) < 1, a relation which can easily be seen geometrically (i.e. ) 

In the random case, where the aj are vector valued random variables on with Lebesgue measure, and i'-"-' is 
a finite dimensional nx n approximation of L and where logdet(L("-') converges to to the Fuglede-Kadison 
determinant logdet(L), we obtain from Hadamard's determinant inequality and Birkhoff's ergodic theorem 

logdet(L) < / log\\a{x,y)\\ dxdy . 

In the case of random Jacobi operators L = t + t* + b{x) = r + t* + Acos(a;) + E, where a{x,y) = 
(...,0,0,1,6(2;), 1,0,...), this gives 

logdet(L) < / log ^2 + (-B + Acos(x))2 dx = log(A/2) + / log v/87A2T4(£7aTcms(^ , 

= log(A/2) + 0(A-i) . 




Improved smaller constants C„(A) < C„_i(A) with Ci(A) ~ C{X) could be obtained by estimating 

/ \og\\A"ix,y)\\dxdy 

from above, where n is small. But already for n = 3, explicit bounds become more difficult due to the 
transcendental nature of the integrals. For n — 2, one would have to estimate the integral 

f , f E ~ Acos(w) — 1 \ / E — Xcos(x) — 1 \ i, , , 

y,.i°gH I )[ 10 jii^^^^ 

from above. With a numerical evaluation of this integral, we get an upper bound for the Lyapunov exponent. 
The values of C„(A) could be estimated computer assisted. We expect that like this, one should reach A„ close 
to 2 for large n. Maybe one could even get A„ ^ for n ^ 00, if also the lower bound of the determinant 
J, I logdetL-ri^Tg^^ could be improved further with the Jensen formula. 



7 Uniform hyperbolicity: a result of Ruelle 

In this section we look at a class of uniformly hyperbolic cocycles and apply a result of Ruelle jl24| to conclude 
that the lower bound on the Lyapunov exponent depends in our case smoothly on parameters. Furthermore, 
the result implies that the neighborhood around such a cocycle, for which one has positive Lyapunov exponents, 
becomes large if the Lyapunov exponent is large. 

Consider a (d + 1) x (c? + 1)- matrix-valued map 

where a{x, y) £ C, d{x, y) is a dx d matrix and b{x, y) and c(x, y)^ are vectors in C*. All entries are assumed to 
be bounded and measurable functions. Together with T € 3^, it defines a cocycle (x, y) i~+ A" (a;, y). We look at 

cocycles in a neighborhood of [ ^ ) , which has one Lyapunov exponent log(A/2) and {d — 1) Lyapunov 
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exponents — oo. 

The maximal Lyapunov exponent in such a uniformly hyperbohc situation (which generalizes the logarithm of 
the Perron- Frobenius eigenvalue of a positive matrix) depends smoothly on parameters jl24| . We formulate this 
in a special case: 

Lemma 7.1 (Ruelle) The maximal Lyapunov exponent A i-^ m(^) o/ is Frechet dijferentiable and realanalytic 
in a neighborhood of Aq — ( J? ) G L°°(T^, + 1,C)), where M{d+ 1,C) is the vector space of all 



X ^ 

complex {d + 1) X (d + 1) matrices. If w ^-^ Aw be a complex parameterization of the cocycle such that for all 
{\w\ < r}, the inequalities 

\aw{x,y)\, ||6u.(a;,?/)||, ||c^,(2:,?/)||, < e 

hold. Then, for fixed A and small enough \e\, or for fixed e > and large enough X, the Lyapunov exponent 
w G {|w| < r} 1-^ fj.{A{w)) is a harmonic function in w. 

Proof. If |o(a;,?;)|, ||6(x,?/)||, ||c(a;,y)||, ||(i(a;,y)|| < e and e is small enough, the cocycle maps the constant cone 
bundle 

ix,y) ^ Cix,y) - C - {c-{l,v) I M < l,cG C } 

strictly into itself. 

(One could write the complex matrix cocycle as a real matrix cocycle on R,2((i+i) ^^j. -^j-^j^j-^ ^y^q cone 

C = {c-(1,?7) I IIt^II <l,7yGR2'*+\cGR} 

in R^'^+^ is strictly invariant. This is equivalent to the fact that one can find a basis and S > such that 
|[A(a;)]ij| > S. By ( 124 4.7), one can deal directly with the complex case. ) 

|124j is written in the case of continuous cone bundles and homeomorphisms, but remark 4.9 in |124j extends 
this to the measure theoretic case. 
From Proposition 4.8 in |124| we obtain especially: 

Corollary 7.2 For small r, the Lyapunov exponent (w, z) ^ ^{G{w, z)) is pluriharmonic (harmonic both in w 
and z) on {(w, z) | < |it;| < r, < [zj < r }. 

Corollary 7.3 Given R > and e > 0. Let Bji{A) be the ball of radius R around A. For large enough A, the 
ball Bfj{G{w, z)) is contained in the open subset of uniformly hyperbolic cocycles and |/i(-B) — /j(^)| < e for all 
B^Br(A). 

Proof. The Frechet derivatives of the function A ^ /^(^) can be computed (see formula (3.6) in |124| ) and 
satisfy D^iJ,{A) < C^/A", where C is a constant independent of A. 

Actually, one only needs the first derivative: Bloch's theorem which says that if a function / is analytic in a 
disc Bfj of radius R, then f{Bfi) contains a disc of radius \Df(a)\Rl, where / is a constant not depending on 



m 



anything (see |3] p. 39). 



8 Determinants of random difference operators 

In this section we look at determinants for classes of random operators. We will see that for a class of ran- 
dom operators, which do not need to be invertible, the determinant is nonnegative and satisfies det[LK) — 
det(L)det(ii'). This product formula of Fuglede and Kadison UB" for operators in a type IIi von Neumann 
algebra usually assumes the operators L and K to be invertible. In the case of random, finite order difference 
operators, we can approximate the operators by periodic operators, which have the spectrum on a union of 
one-dimensional curves. By changing the energy we can approximate the operators by invertible operators for 
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which the formula holds. The formula holds then for the random difference operators also. 

A measure dk in the complex plane is called the density of states of a (not necessarily selfadjoint) random 
bounded operator {x,y) G i-^- L(x,y) = X]j=-p (^7 y)''"'' /^(Z,C^) with N x N matrices a'^^^ and 
TUn — Un+ii if for every polynomial /, one has 

/ Trace([/(L(a;,y))]oo) dxdy = / I{E) dk{E) 

where Trace is the usual trace for iV x iV matrices normalized in such a way that Tracc(l) = 1. 

For operators L with a density of states dk one has a trace 

tr(L) = j E dk{E) = J Trace[L(x, y)]oo dx dy . 

If the Thoulcss formula holds and the Lyapunov exponents of the transfer cocycles are all finite, then the mea- 
sure dk has finite logarithmic energy and the potential E 1-^ tr(log \L — E\) = log \E — E'\ dk{E') is defined 
and finite for all E. 

In general, if L has a density of states dk, a determinant det(L) of L can be defined by 

logdet(L) = tr(log(|L|)). 
If the value J^-, log \ E\ dk{E) of dk at the point E = should be —00, one would put det|L| = 0. 

(The extension of the determinant to noninvertible operators, which we look at, is what Fuglede and Kadison 
call an 'analytic extension of the determinant' to singular operators. This is in contrast to the 'algebraic exten- 
sion' for which the determinant is put always to zero if the operator is noninvertible). 

Given T,S ey and a^^^ , 6^^) S (T^ , dxdy) . Two not necessarily selfadjoint elements L = a^^^ i.^, y)T^ , K = 
b^^\x, y)a^ in the crossed products define random operators 

p p 
{L{x,y)u)n= ^ a^^\x,y)un+j, {K{x,y)u)n= ^ b^^\x,y)un+j 
j=-p j=-p 

with a'^\x,y) = a(J^T''{x,y)),b'^\x,y) = b^J^S^ {x , y)) . 

Lemma 8.1 The operators L, K, LK have all density 0/ states. The Thouless formula holds for L, K and LK 
and det{LK) = det(i)det(i^). 

Proof, (i) The operators L, K can be embedded in a single von Neumann algebra without changing the deter- 
minant. 

Proof. L is in the crossed product Zt of L°°{T^,dxdy) with the Z-action ar ■ f f{T) and K is in 
the crossed product Z5 of L°°{T^,dxdy) constructed with the Z-action as '■ f ^ f{S). The operator LK 
is in the crossed product ZtxS of L°°{T^ x T^,dxdy x dxdy), constructed with the Z-action arxs, where 
T X S{xi,yi,X2,y2) — {T{xi,yi), S{x2,y2)) on x T^. One can embed the two crossed products Zt and Zs 
in Ztxs- The density of states of LK exists in the same way as it exists for L and K. Moreover, the density 
of states of L (rsp. K) in ZtxS is the same as the density of states of L in Zt (rsp. in Zs). 

(ii) For the verification of the Thouless formula, see Section 15). 

(iii) We can assume without loss of generality that L and K are invertible. 
Proof. By ergodic decomposition (i.e. [2^1)) one can assume T, 5* to be ergodic. 

Approximate T,S G 3^ by periodic transformations Tk using Rholin's lemma |53| : the measure of {{x,y) € 
T^, T{x,y) 7^ Tk{x,y)} is then smaller than 1/fc and T^ = Id. By ergodic decomposition, the determinant 
is in such a case an integral of determinants of type Hi factors, in which case Lk,Kk as well as the product 
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LkKk are periodic operators which have their spectrum contained in a finite union of real smooth curves in 
the complex plane. In this case, there are arbitrarily small complex numbers Ek such that Lk — Ek, Kk — Ek 
are invertible allowing to apply the determinant formula. Now, E i— > det(Lfc — E),E i-^ det{Kk — E) as well 
a.s E 1-^ det((Lfe — E){Kk — E)) are continuous. Lk — Ek,Kk — Ek are invertible and det{Lk — Ek) det(L), 
det{Kk - Ek) det(A'), and det((ifc - Ek)iKk - Ek)) ~^ det((i - E){K - E)) for n ^ oo. 

(iv) The product formula. 

Proof. We can now invoke the determinant theory for finite type von Neumann algebras |2H1 (see |27| chapter 
I) which tells that det{LK) — det{L)det{K) for invertible operators L,K in a. type IIi von Neumann algebra. 



Corollary 8.2 If L^^"^ = a^*V+ (c^'V)* +b^^^ are random Jacobi matrices over dynamical systems T^"^^ : 
then 

det(L(i)x(2)) ^ det(x(i)) det(L(2)) . 



Proof. Apply Lemma H8.1|l . 



Whatever Lyapunov exponent estimates for an analytic cocycle can be proven for irrational rotations a; i— *■ x + a 
as in |56[I133| . they can up to some correction C(A) also be proven for this cocycle, where the underlying system 
is an arbitrary measure preserving transformation T ^ y. For example, for any realanalytic nonconstant /, the 
Lyapunov exponent of AT^xf should be positive for large enough A. 



9 The Jensen formula in a sector 



In this section we prove an integrated form of the Cauchy-Riemann differential equations. Unlike the classical 
Jensen formula which holds in an annulus (for a short proof see 43 ), we formulate it in a sector. This allows 
to apply it in piecewise analytic situations. While we don't use this formula explicitly, we used it to explain 
some heuristics in the introduction. We furthermore think it could be useful for further improvements on the 
Lyapunov estimates. 

Consider a sector ^ = {z e C | |2:| € 
[r, i?], arg(z) g [a, fi] } and denote by "ft = {\z\ — 
t} r\ A a, circular arc in A. Let g be a bounded 
continuous, complex- valued function on A which 
is analytic and nonzero on ^ \ 5, where 5 is a 
finite set of points in A. Define for t G [r, i?] 

Arg^^(ff) = Re— / ^(z) dz . 

(It is defined except possibly at finitely many t, 
where 74 fl 5 is nonempty.) 

Lemma 9.1 (Jensen formula in a sector) If g is analytic in A and nonzero in A\S, where S is a finite 
set of points, then 

dt 

/ log |g(i?e2-")| dx = / log \g{re'^n\ dx + / --Arg^^(g) dt . 
Proof, (i) Assume first that g is analytic in A and that it has no roots in A. Clearly, 




R 



log \giRe^^n\ dx = log \gire^^")\ '^^ + ^ ^ l5(te'"")| dt dx 



Because g is analytic, the Cauchy-Riemann equations -^g{z) 
{d/dz)g{z), we obtain 

g'{te^^^^) 2Trix 



d_ 

dt 



log|5(te^"-)| =Re 



g(te2^^^) 



hold. With z 
' g'{z) z' 



t exjp{2TTix) , g'{z) 



Re 



9i^) t 
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We have therefore, using 2nidx = dz/z, 

/ log|5(i?e2"-)|dx- / log|5(re2-")|dx = Re—/ - dt 

Ja Ja 27rj J g[z) 



z) t 



dt 
T 



(ii) In general, we partition A = [f , -R] x into smaU sectors ^^fc = [rj,rj_|_i] x [Q;fe,Q;j.+i] and apply the 

formula to small sectors Ajk which do not intersect 5. The assumption on S and g assures that for all the 



integrals, the contribution of the sectors which intersect S can be neglected in the limit n ^ oo. 



Remarks. 

1) The Jensen formula in the annulus \l'6'6l llUj follows: if A is the annulus and if g is analytic in {arg(z) G 
(ai, tti+i), \z\ £ (r, R) } and if g{z) = YVjLii^ ~ o,j)h{z) such that h{z) ^ everywhere on A, then 

/ log|ff(i?e2-")|dx- / \og\gire'^m dx = Arg^^ (/i) + V log( — ) . 

Jo Jo Jr I PjI 

2) The Jensen formula for subharmonic functions in the annulus A follows: for subharmonic functions z <—>■ 
log 13(2)1 = / log \z — z'\ dk{z') in the disc {\z\ < R}, one has 

r log dx = r log \g{re'^n\ dx + f log{^) dk{z') . 

Jo Jo JAr,R 1^1 

3) Lemma (|9.1|) holds for sets S which are bigger than just a finite point set. For example, S can be a subset 
of a finite union of smooth curves. 



10 Toral homeomorphisms and Lax approximation 

In this section we review an approximation result for homeomorphisms on the torus. It allows to approximate 
the group of measure preserving homeomorphisms by a finite group of piecewise linear transformations. In- 
vestigating the Lyapunov map T 1-^ A*(^T,Asin) on finite groups of measure preserving transformations was 
important for us because the evaluation of the Lyapunov exponent is a finite dimensional integration and quite 
reliable. 

On the group y of all measurable, invertible transformations on the torus T*^ which preserve the Lebesgue 
measure, one has the metric 

p(Ti,T2) = |dist(ri(x),T2(a:))|oo , 
where dist is the geodesic distance on the flat torus and where | • |oo is the L°°-norm. 

A cube exchange transformation on T'' is a periodic, piecewise afhne measure-preserving transformation which 
permutes rigidly all the cubes Yli^iih/n, {ki + l)/n], where fc^ G {0, . . . , n — 1}. Every point in T'' is periodic. 
Such a transformation is determined by a permutation of the set {1, . . . ,n}''. If the permutation is cyclic, the 
exchange transformation is called cyclic. A theorem of Lax states that every T £ X can be approximated 
in the metric p by cube exchange transformations. Actually, Tk can be chosen to be cyclic 

Lemma 10.1 (Lax approximation) For every e > and every homeomorphism T £ X , there exists a cyclic 
cube exchange transformation which satisfies p{Tk,T) < e. 

Lax's 'book proof of this result uses Hall's marriage theorem in graph theory (for a 'book proof of the later 
theorem, see [J). 

Remarks. 

1) Periodic approximations of symplectic maps work surprisingly well for relatively small n (see |12U) 1. On 
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the Pesin region this can be explained in part by the shadowing property |65j . The approximation by cyclic 
transformations make however long time stability questions look different |51) . 

2) It has been measured that certain Lax approximations T„ of the Standard map have longest orbits of size 
> Sn"^ with (5 > 

3) If r G y^, the density of states dk{x, y) of a periodic operator L{x^ y) is contained in a finite union of curves. 
The density of states of L is the measure /q^^'' Jq^'^ dk{x, y) dxdy which has positive area in the complex plane. 
For a cyclic Tk, the Lyapunov exponent is computed as 

.i/k ni/k ^ 

li{ATj = k-^k^ / log X^^^iA'' {x,y)) dxdy , 

Jo Jo 

where Xmax (A) is the maximal eigenvalue of a matrix A. The factor k^ comes because we integrate over a region 
with Lebesgue measure 

11 Thouless formula for nonselfadjoint operators 

In this section we look at the Thouless formula for random operators which can have spectrum with positive 
area in the complex plane. Nonselfadjoint operators appear naturally, when a problem is parameterized by an 
external complex parameter like for operators . The proof of the Thouless formula which we are looking at in 
this section relies on the fact that for a general parameterization, there is an abstract Thouless formula which 
relies on the Riesz decomposition theorem in potential theory. 

Consider an analytically parameterized, not necessarily selfadjoint, random Jacobi matrix 

(6) 

with a complex potential parameterized by a complex parameter z. Examples are the w- parameterization 

fw{x) = (w^"'" exp(ia;) + w exp(— i.-E))/2 (7) 

for C \ {0} and for fixed w ^ Q, the energy parameterization 

fEA^.x)^ f^{x)-E . (8) 

Proposition 11.1 (Thouless formula for nonselfadjoint operators) In the energy parameterization case 
f^). the Riesz measure of the Lyapunov exponent ^(At.e) of the transfer cocycle At.e is the density of states 
of the operator 0). 

Proof. For any general analytic parameterization z i— s- L{z), the Lyapunov exponent ^i{Az) of L{z) is a a 
subharmonic function of z. The Riesz decomposition theorem (see i.e. |54l \121\ ) implies the abstract Thouless 
formula = / log \z — z'\ dk(z') + h(z), where z i—f h{z) is harmonic. 

The name 'abstract Thouless formula' is from |56| . 

If z i-^ ^{z) is harmonic for large \z\ and grows like log|az| for \z\ i—f oo then ij,{z) = log(a) + Jlogjz — z'| dk[z'). 

The Riesz measure dk of the subharmonic function z i-^ pl{z) satisfies dk = in the sense of distributions, 
where A is the Laplacian on distributions. The measure dk coincides in the periodic case with the density 
of states because in this case the Lyapunov exponent is zero exactly on the spectrum which implies that the 
density of states is the equilibrium measure on the spectrum. 

Consider the Riesz measures dfc„(z) — A/i„(z) of the subharmonic functions 

z ^ ^J.n{z) ^ n^^ / \og\{A'lj.{x,y))ii\dxdy 

JT2 
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and let dk{z) — A/i(z) be the Riesz measure of the Lyapunov exponent z i— > ^(^Az^t)- Because ^n{Az^T){x,y) 
converges by the multiphcative ergodic theorem to ^{Az^T){x,y) for almost all {x,y), the convergence holds 
also in the sense of distributions. Then also the Riesz measures c?fc„ = A/i„ converge to A/z as distributions 
for n oo. But then, because smooth functions are dense in all continuous functions, dkn = A/z„ converges 
weakly as measures to dk = A/i. 

In the energy parameterization case (jSJ , the density of states of L is defined as the measure dk in the parameter 
plane C satisfying Jrj,2[g{L{x,y))]oo dxdy = j g{E') dk{E') for every continuous function g. 

For E outside the spectrum of L, the sequence Jq;{E' ~ E)^^ dkn{x, y, E') converges for n ^ cxd uniformly to 
J^{E' - E)-^ dk{x,y,E'). Because 

tr((L(")(x,2/)-i?)-i)- / iE'-E)-Uk^ix,y,E'), 

Jc 

where the Riesz measure dkn{x,y) = A/i„(x, y) is a finite point measure located on the point spectrum of 
defined in an earlier section, the formula holds also after integration over and dkn is the density 
of states of the random operator (x, y) i-^ L^'^\x, y). 

The functions gj(E) — [E — Zj)~^ , where {zj}jgN is a dense set in the complement of the spectrum of L, span 
a dense set in C(C). For such functions, we have 

/ [9j{L^''\x,y))]m dxdy ^ / [gj{L{x.y))\ao dxdy 

by the Avron-Simon lemma. 

We know therefore /.pa [^(-^(s;, ?/))]oo dxdy — J g{E') dk{E') for all g £ C(C) so that dk is the density of states. 




Examples. 

1) If w 1-^ B{w) = 'w''A{w) is analytic in C, where A{w) is the transfer cocycle of L^, then the Lyapunov 
exponent of B satisfies 

fi{B{w))=C+ I log\w -w'\ dk{w') 
Jc 

where dk has total mass + 1) and where C is a constant. 

2) For periodic operators, where the support of the Riesz measure is on spectral bands, the Lyapunov exponent 
oi w 1-^ A{w) is = Re[i arccos(Tr(A")(w)/2)]. With the parameterization 10, we have 

= -log|w| + / \og\w - w'\ dk{w') +log\\/2\ 
Jc 

because w i-^ Iji{A{w)) grows like log|wA/2| for \w\ — > oo and dk{C) = 2. The Lyapunov exponent of B{w) = 
wA(w) which is harmonic near satisfies then 

^{B{w))= I \og\w - w'\dk{w') + \og\X/2\ . (9) 
Jc 

We expect the function w i— s- fi{B{'w)), which is harmonic near w — to extend often to a harmonic function in 
the entire complex plane. 

Remark. By the Hadamard three circle theorem the function fj.maxi'r) = max.\^\^r Im(^)I is a increasing 
function in r which is convex in log(r). Aiming to get rid of the C(E, A) term in Conjecture Hl.l|l . it would be 
good to know M — {\w\ — 1 | — fJ-maxi\w\) } which is contained in the subset of the unit circle on which 

the potential is > u(0). In the case of the Standard map, a candidate for a maximum on |w| = 1 is w = 1 
because the pointwise Lyapunov exponent ii{x,y,w) satisfies iJ,{x,y,w) — ij{—x,—y,w) so that fJ,{w) = /i(w). 
See 18) in the discussion section. 
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12 Log-Holder continuity and capacity 



In this section we consider a potential theoretical property of the Riesz measure of the Lyapunov exponent. 
This applies to the Riesz measure of the parameterization w i— > or, for fixed w to the energy parmeterisation 
E 1-^- L^u — E, where the Riesz measure dk is the density of states. 

As in the selfadjoint case, the fact that the Lyapunov exponent is nonnegative implies mild regularity for the 
measure dk. Log-Holder continuity is a property of measures with finite potential theoretical energy. It implies 
that such a measure is absolutely continuous with respect to the zero-dimensional Hausdorff measure |35) . 

Lemma 12.1 (Log-Holder continuity) Letw > L^, be an analytic parameterization of a random Jacobi op- 
erator. The Riesz measure dk of the nonnegative subharmonic function w '—^ fi{AT{w)) is log-Holder continuous: 
there exists a constant C such that 

dk{B)<C -{logil /\B\)r' 
for every ball B with diameter \B\. The constant C only depends on L and \E\. 

Proof. The proof is the same as in the real case |22 El- assume dk is contained in a ball of radius c and 
assume that B is any ball of radius r < 1 around z S C. With A = {z' \ \z' ~ z\>\ }, we have 

0< / \og\z - z'\dk{z') <\og{l + c+\z\) . 

J A 

Therefore, using in the first inequality of the following identities also that |z — z'| < 1 on and log \ z — 
z'\ dk{z') > 0, we have 

0> / \og\z- z'\dk{z')= [ \og\z- z'\dk{z')~ j \og\z - z'\ dk{z') > Q ~ \og{l + c + \z\) . 

J A" JC J A 

For zi ^ z with \zi — z\ — r, and because \z — z'\ < \z — zi\ for z' e 

\og\z - zi\ dk{B) > I \og\z- z'\dk{z')>^\og{l + c+\z\) 

J A' 



SO that < log(l -t- c\z\) -f log \z — zi\ dk{B) which is the claim with C — log(l 4- c\z\). 



Example. Consider the case, when T = Tfe is a Lax cube exchange transformation on the torus and where the 
operator 

Lw{x, y) = T + T* + X/2{w^^Zj wzj^) 
with Zj — Zj{x, y) — exp(ia;j), {xj, yj) — T^ix, y) is a periodic, nonself adjoint Jacobi matrix. 



/ 



We are interested in the potential theory 
of the w-spectrum 

a(L(x,y)) = {«;|tr(A^(a;,y))e [-2,2]}, 



which supports the 'w-density of states' 
I dk{x,y) of L{x,y). In the picture to the 

^ left, we see an example of the 'w-spectrum' 

\ J a{L{x,y)) for A = 2.1 where k = 7. One 

^ can see six bands located on the unit circle 

and eight bands away from the unit circle. 

The 'w-density of states' of L is then dk ~ dk(x, y) and its support, the 'w-spectrum' of the random 
operator L is the union IJ^ yex^ '^{L{x,y)) and has in general positive area in the complex plane. 
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The potential theoretical energy 

lidkix, v)) = - ( ( log \z - z'\ dk{x, y)iz') dk{x, v){z) 
Jc Jc 

is I{dk{x,y)) — log(A/2) dk{x,y){C) — 21og(A/2) because dk{x,y){C) ~ 2 and the Thouless formula assures 
that the Lyapunov exponent ^{x,y){z) is related to the potential u{x,y){z) = J(-,log|z — z'\ dk{x,y){z') of 
dk{x,y) by ^{x,y){z) — log(A/2) + u(x,y)(z) — log{'w) and the Lyapunov exponent is zero on the spectrum 
a{L{x,y)). The measure dk{x,y)/2 is the equilibrium measure on the set a{L(x,y)). It has the potential 

theoretical energy log(A/2)/2 and the capacity C{a{L{x,y))) is J j-. 



13 Harmonic continuation of some singular integrals 

In the case of the Standard map, the function fi{w) ~ log(ti; det(Lm)) with given by is for small \w\ 
the real part of an analytic function. By the Cauchy-Riemann differential equations, the arg(i/;)-dependence of 
the function w i—s- exp(/i(u))) is related to the change of the harmonic conjugate along the radial direction. The 
Jensen formula extends this to the t«-spectrum of L^, where /i is no more analytic. There are cases (and we 
expect it to be the case in the Standard map), in which the harmonic conjugate can be extended as a harmonic 
function on a neighborhood of the unit disc. Estimating this harmonic conjugate with complex analytic meth- 
ods could lead to possibly better lower bounds of the Lyapunov exponent. Alternatively, one could look at the 
argument change (which differs from the harmonic conjugation) and use the Jensen formula directly to estimate 
the oscillation of the Lyapunov exponent w i— > /i(yl(w)) for w the unit circle. In this section we formulate 
the problem of computing the harmonic continuation rsp. the argument change along radial lines, when the 
w-spectrum is located on the unit circle. In the next section, we will see that assuming dk to be supported on 
T is no big loss of generality for estimation purposes because projecting the spectrum onto the unit circle does 
not change much the potential (which is the Lyapunov exponent.) 

A continuous, real-valued periodic function x ^ fpi^) on the circle defines a complex-valued function x i—f 
a{x) = exp(i(a; -I- (jji^))) ^-nd a measure dk = h{x)dx on the circle T by pushing forward the Lebesgue measure 
dk = a*dx. Denote by a„ = J exp(—inx) h{x)dx the Fourier coefficients of this measure. The function 
/(z) — J!p dk{x')/{z — exp{ix')) — Jr^ dx/ (z — a{x)) defines two analytic functions g±{z), where g+ is analytic 
on |z| < 1 and 5-(z) is analytic in |z| > 1. These functions are defined by their Taylor expansions at z+ = 
and z_ = oo: 



h+iz) = -y^a-i(a;)^— ^da: = -X:z"a_(„+i), 



„ ^ OO —OO 

h^{z) = / ^-i^^^^da: = ^z-("+i)a„= ^ z"a_(„+i) 



n=0 

The function h = — = X^nez •^^"^-(n+i) is like this represented by the boundary value of two analytic 
functions h±. One has for example that if h is in L^, then h± are in the Hardy space Rp. The function 
h is realanalytic if and only if h is analytic in some strip around the unit circle. In this case, we can write 
/(z) = h+{z) = h{z) + h_{z) for {1 - e < |z| < 1} and /(z) = /i_(z) = -h{z) + h+{z) for {1 < |z| < 1 + e}. 
The functions h± can then be analytically continued to a larger domain. If dk extends to an entire function, 
then both h± extend to entire functions. 

The integral H{z) ~ Jr^ log(z — a{x)) dx with H(0) — can be expanded in the same way at z = and oo and 

H^H_-H+^J2 ^"+'«-(n+i)/(" + 1) = E ^"«-"/" ■ 
nez nez 

The imaginary part of on |z| = 1 is the average argument change from to oo of the function a z — a(x) 
along the line going through z = exp(jQ;). 
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Near :— 0, the function Im(i7_(2;)) is the argument change of exp(77_(z)) from to z. Near Zj^ :— 0, the 
function Im(_ff_|_(z)) is the argument change of exp(iJ+(z)) from oo to z. 

Because of the Jensen formula (I9.1II . we are interested in the argument change a(2;) from to z along a straight 
line of the function Jrplog(2; — a(x)) dx. For each point a(x), the argument change a(z) from to z is a well 
defined angle. For \z\ < 1, we have a{z) — lm{H_{z)). For \z\ > 1 we can write a{z) — Im(iJ+(z)) + /3(z), 
where /3{z) = (3{z/\z\) is the argument change along the straight Hne from to cxd which passes through z. Now 

/3(z) = lm{H{z/\z\)) so that 

a{z) = lmiH+{z)) + H{z/\z\)) , \z\ > 1 
a{z) = Im(i/_(z)) , |z| < 1 . 



14 Lipshitz approximation of the Lyapunov exponent 

In the theory of Schrodinger operators one uses the fact that the Lyapunov exponent ij-{Ae) of the transfer 
cocycle Ae of L is the Hilbert transform of the integrated density of states (or rotation number [1231 \2b\ ) 
p{E) = dk{—oo,E), where dk is the density of states of L. For the spectral problem w i-^ L^, where the 
w-spectrum is contained in an annulus Ar = {\z\ G (r, r^^) },r < 1 containing the unit circle, one can ask for 
the regularity of the map a i— *■ //(L(exp(ia)). The situation for the unit circle in the last section is relevant to 
this problem. The idea is to project the measure dk onto the unit circle with the map 7t{z) = arg(z). While 
this does not change the Lyapunov exponent much, one can use the classical properties for the new Lyapunov 
exponent which is the potential of a measure on the unit circle. A Lipshitz continuity of the 'integrated density 
of states' a i— > dk{{z \ arg(z) g (0, a)}) would imply that a i-^ ii{a) — yLt(L(exp(ia))) would be close to a 
Lipshitz continuous map. 

Let dk be the Riesz measure of the Lyapunov exponent of w i— > L^j so that = /^^ log |z — z'| dk{z'). Let 
■n*dk be the measure on the unit circle obtained as the push-forward of dk under the projection 7r(r exp(i(/))) — 

The fact that the Hilbert transform preserves Lipshitz continuity can be derived from the Jensen formula in 
a sector (rsp. the Cauchy-Riemann differential equations). Let H be the function as in the last section 

with the measure 'K*dk on the unit circle. The statement is that if the radial argument change lm{H{z)) of the 
function z ^ g{z) = exp(iJ(z)) is Lipshitz continuous on the unit circle and the Lipshitz constant is C, then 
the function Re(i?(z)) = / log \z — z'\ ■n*dk{z') is Lipshitz continuous on the unit circle with Lipshitz constant 
C. 

Lemma 14.1 // dk is supported on Ar — {\z\ G (r, r^^)} and dk is the Riesz measure of a Lyapunov exponent, 
then 

\og\z - z'\ dk{z') > I \og\z~ z'\'K*{dk){z')-d{r) . 
c Jc 

where dir) — > for r — 1 . 

Proof. Given a point zq with \zq\ = 1. The set Y{zQ,r) = {w £ Ar \ \w — zq\ > \Tr{w) — zo| } is an open set. 
We have F(zo,r) C F(zo,r') for r < r' . Write 



\og\z ~ z'\ dk{z') ^ \og\z- z'\dk{z')+ \og\z - z'\ dk{z') 

C JY{za,r) JAr--Y(zo,r) 

We have /x-y(zo,r) 1^ ~ ^'1 *(^') - lAr-Y{zo,r) - z'\ 7r*dfc(z') and 



\og\z - z'\ dk{z') > / log|z- z'l 7r**(^') -^^C'-) 

Y(zo,r} JY(zo,r) 



where d{r) 0. 

A numerical bound on the Lipshitz constant of a i-^ n* {dk){f3) d(3 = dk{{z \ arg(z) G (0,0;)}), where dk 
is the Riesz measure of the Lyapunov exponent w ^{AE,\f,Tiw)) would allow to estimate the Lyapunov 
exponent of the transfer cocycle ^_b,a/,t of the operator Q). 
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15 An estimate for harmonic maps 

We have numerical evidence that in the case of the Standard map the Lyapunov exponent fi{w) is the sum of 
a positive subharmonic function and a harmonic function obtained by harmonic continuation of Re J log{'w — 
w') dk{w') + log |w| from a neighborhood of {w = 0}. If this should turn out to be true, estimates for harmonic 
maps could be used to give alternative and maybe better estimates for the Lyapunov exponent. 

For a general harmonic function h, one can estimate inf i^i^^ h{z) in terms of C(r) — supi^i^^ h{z). We present 
two ways. The first way uses Harnack's inequality. A second approach uses the distortion theorem for univalent 
functions. 

Lemma 15.1 (Lower bound for harmonic maps) Let h be a harmonic map in a disc D(r') = {|z| < r'} 
with r' > r,r' > 1, such that h{0) — and h{z) < C(r) for all z G D(r). Then a) 

h{z) > ~C{r)^— . 



for \z\ — 1 and b) 

(1+r-^) 
Mz)>-C(r)(l + ^^-^) 

for \z\ = 1. 

Proof, a) v{z) = C(r) — h(z) is nonnegative in the disc D(r) = {\z\ < r} and v{0) — C{r). Harnack's inequality 
gives 

C(r) - h(z) = v(z) < C(r)^-^ 

r — 1 

so that h{z) - C{r) > -C(r)^ and h(z) > C{r){l - ^) - -C(r)^. 

b) The function C(r) — h{z) is positive in D(r). Let H(z) be the analytic function in D(r) which has h as its 
real part and which vanishes at 0. By the Noshiro-Warschawski criterium |33| . the analytic function 



F{z)= r{C{r)-H{z)) 
Jo 



dz 



is univalent in D(r) because it has a derivative which has a positive real part in D(r). The function F{z) = 
F{rz)/{rC{r)) is in the class S of functions which are analytic and univalent in the unit disc and satisfy 
F'(0) = 1. By the distortion theorem for maps in the class S, we have 

in.)!**'""' 



(1 - R)3 

for \z\ ^ R < 1. Assuming r > 1 and applying this for R = r^^ , we get for \z\ ~ R 

' ^ ' C(r) ' ' C{r) ' - (l-r-i)3 " 

In other words, if \z\ = 1, then 



' C(r) ' - (1 -r-i)3 
which leads to < C(r)(l + ^^/^Jlip ) and implies the claim. 
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16 Remarks 



1) Measure preserving maps. 





S,T,,4cos) - log III 


0.2 




^ fTn 




/ -0.2 





All cyclic permutations of {1, • ■ • , 9} . They 

are ordered such that the excess to 

log |A/2|, shown here for A = 4, increases. 



The random variables Xn where T"{x,y) — 
{xn, Hn) have as their law the Lebesgue measure. 
One could enlarge the class of transformations 
even so it is not much gain of generality |87| . 
The conjectured lower bound log(A/2) can be 
expected to hold in general only for homeomor- 
phisms T e A". The reason is that for some 
piecewise periodic measure preserving transfor- 
mations : {x, y) ^ (x -\- Tr{x)/n mod 1, y) with 
a permutation tt on {l,2,...,n}, the value of 
m(^t,Acos) can be slightly smaller than log(A/2). 
The Figure to the left shows the distribution 
of the random variable tt i— > /^(vIt^.Acos) on the 
group 3^9 C 3^ of all permutations of the 9 annuli 
{(i/9, 0' + l)/9] X T C T'^}%i in the case A = 4.0. 



It would be interesting to know inirp^yfi^Ax^xcos) for given A > and whether the infimum is attained. By 
upper continuity of the Lyapunov exponent the infimum taken on the discrete set |J„ 3^„ C 3^ is the same as 
the infimum taken on X. 



2) Hamiltonian flows. 

Like for a large class of monotone twist maps, any Standard map is the time-one map of a Hamiltonian 
differential equation with time-dependent periodic Hamiltonians This is also true in higher dimen- 

sional cases (111 i49] . If positive Kolmogorov-Sinai entropy is stable under perturbations of the map, there 
would be many real-analytic, time-dependent periodic potentials V for which the time-dependent Hamilto- 
nian system x = V{t, x) has positive Kolmogorov-Sinai entropy. To see this, we note that the potential 
V{t,x) = 6ti£z^smm{x) (which is a distribution in t) has as a Poincare map the Standard map TAsin™- The 
smoothed-out potentials Ve{t,x) = (peit) * V(t,x), where (pe is a mollifier function, has a Poincare map which is 
realanalytic in (x, y) and which is arbitrarily close to Txsinm in a Banach space of realanalytic maps. For fixed A, 
the Poincare map has now positive metric entropy for small e and large m because the stability result as stated 
in Conjecture extends to general realanalytic perturbations of the map (and not only to perturbations of 
the function /). It would be nice to see entropy estimates for Hamiltonian flows on three dimensional energy 
surfaces of systems x = V{x), where V{x) is a potential on or to geodesic flows on the torus or to real- 
analytic, strictly convex Birkhoff billiards for which we have no example, for which positive Kolmogorov-Sinai 
entropy has been proven. 

Even more off-limit seem problems in celestial mechanics like the Stormer problem of particles trapped in the 
van Allen belts of the Earth's magnetic dipole field Jl] (for which Pesin theory in the form jSH] can be shown 
to apply) or particular Newtonian three body problems. 

3) Smaller A. 

Estimating the entropy in the case of small A is also open for the Chirikov Standard map. Experiments indicate 
that /i(A) > for aU A > 0. 

Numerically, we find /i(A) > log(A/2) for all A and as a "rule of thumb" that the averaged Lyapunov exponent 
^(A) satisfies almost but not exactly the Aubry duality /x(A) ~ log(A/2) -I- /.t(4/A). 
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/i(^Acos)-l0g+|f| 



0.4 



4/A 10 8 6 



6 8 10 A 



How close are we to Aubry duality? 

For each A we averaged 70x70=4900 different 

orbits of length 10*" on the phase space. 



The duality holds in the case of the dynamical 
system Tq (the almost Mathieu case). It does 
not hold for T^sin as the Figure to the left shows 
for A e [4/10,2] and A e [2,10]. Indeed, by 
computing moments of the density of states 
of L using a random walk expansion j73j . one 
can check that the discrete random Schrodinger 



operator Lt,Ac 



Acos(a;) 



over 



the 



dynamical system T\ does not have the same 
density of states as the naive "dual operator" 
L - (A/2)A + 2cos(x) = (A/2)Lt,(4/a)cos- (K 
is not excluded however that changing T in L 
could restore a generalized duality.) 



As in the almost Mathieu case, there is a two-dimensional magnetic operator involved. The magnetic field at a 
plaquette Pn,m is exp(i(a;„ — Xn-i) independent of m. In the Standard map case, the magnetic fields depend 
on space and are correlated. See 25) for more on the Aubry duality. 

4) Smooth cocycles. 

Results like in 154, 156_ depend on the dynamical system, e.g. Diophantine conditions for irrational rotations. 
It is known that positive Lyapunov exponents of SL{2, R) cocycles can be destroyed away from Anosov systems 
by measurable perturbations of the cocycle jfi8j , by continuous perturbations 39^ or using the unprovcn " Last 
theorem of Mane" [SS], by C^-perturbations of a non- Anosov map. 

5) Equilibrium measures. 

Via the variational principle (see e.g. |143) ) metric entropy results give lower bounds on the topological entrop. 
with previous known lower bounds |32l I75j . The estimates in 75 are obtained using a Gronwall argument 
which estimates the spectral radius of the operator (where L is defined by a hyperbolic invariant measure 
of the Chirikov map) and which are far from optimal (the bounds in 1751 are less good than the estimates in [221 
but are proven when / is a nonconstant Morse function and not only in the case of the Standard map where 
/ = sin.) The equilibrium measure is a T-invariant measure on which maximizes the entropy and has as 
the metric entropy the topological entropy. In the case of the Standard map, it is not known whether there is 



Lyapunov exaponents have applications in the 
theory of discrete random Schrodinger opera- 
tors El inil- The attribute 'random' is 
used in that theory often in the same way as 
'random variable' is used in probability the- 
ory. Indeed, random operators are classes of 
operator- valued random variables. There is in 
general no randomness assumed in the potential. 
For example, if T is a translation on the torus, 
the random operators are quasiperiodic. Using 
Pastur-Kotani theory. Conjecture ?? enlarges the 
class of ergodic discrete one-dimensional random 
Schrodinger operators L, for which there is no 
absolutely continuous spectrum. 

Conjecture ?? implies with 111 that for a Baire generic ergodic T X, the corresponding ergodic operators 



an absolutely continuous equilibrium measure. 
6) Absolutely continuous spectrum. 
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The Lyapunov exponent of the Chirikov Hamiltonian. 
The graph shows a nonnegative excess to log lA/21 
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have no absolutely continuous spectrum for large A. ^ For non-ergodic T , the absolutely continuous spectrum is 
only absent for a large set of a; g T'' because the Lyapunov exponent can be zero on a set of positive measure. 
Indeed, for invariant measures on KAM tori, the corresponding ergodic almost periodic Schrodinger operator 
can have some absolutely continuous spectrum for small A |78| . ^ The size of the spectrum, where the Lyapunov 
exponent vanishes can be estimated if the orbits do not stay too long at the same x value for a long time |138| . 

7) Discrete Spectrum. 

The question of localization of the operators L\cos |134| stays open. The conjugation of the map to a Bernoulli 
shift on some positive Lebcsguc measure suggests that eigenvalues should exist for similar reasons as in the case 

In the Figure to the left, numerical re- 
sults are shown for the Fourier coefficients 
[In = ((/), exp(— mL)(/)) = J^e^'"^ d^^{9) of the 
spectral measure /j,^ where (j> — (. . . , 0, 1, 0, . . .) G 
P{Z,C). We see the case of the Standard map 
operator and Mathieu operator. With a Wiener 
criterion for the existence of discrete spectrum, 
(where we compute a fast discrete quantum 
evolution using a related operator which has 
the same type of spectral measures [ZSj), we 
find numerical evidence for some point spectrum 
for L{x,y) for a set of {x,y) S of positive 
measure, if A is large. 

In the picture to the left, we confirm numerically 
the theoretically established localization transi- 
tion for the Mathieu operator at A = 2 |6(J[ 1^ . 
In the picture to the right we see an indication 
for some point spectrum in the Chirikov case for 
A > 0. 



Note that in the same way as for the Anderson operator, the operator LT{x,y) has for an Anosov T purely 
singular continuous spectrum for a Baire generic set of {x,y). This can be derived from general principles |13l)| . 
Also, for Baire generic (T, {x, y)) E X x and A > Ao, one knows that Lt, \cos{x, y) has purely singular continu- 
ous spectrum. This means that y) on ^^(Z, C) has Baire generically no bound states and no extended states. 

8) The distribution of the Lyapunov exponents. 

^The space X becomes a complete metric space with the metric p(T, S) + p{T~^ , S~^). 

•^The results in 1781 have not been written down in a final form since both of us got engaged in other projects. 



of independent, identically distributed potentials. 



" ^ Y.k=l lAfc(iAcos)P 
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Wiener test for localisation for AG [0, 10]. To the left, the 
Mathieu case (n=100'000), to the right, the Chiricov case 
(n=10'000 averaged over 50x50=2500 orbits). 
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The distribution functions of the random variable 
for A = 2, 10. For A = 2, m[Xx = 0) > (KAM). 



For every A, the pointwise Lyapunov expo- 
nent {x,y) I— > fi\{x,y) is a random variable 
on which is T-invariant and has a mean 
> log(|A|/2) — C(A). The numerical exper- 
iments indicate that the random variables 
X\{x,y) = pL\{x,y) have the property that the 
normahzations = {Xx - ¥j[Xx]) / ^jYa,r[Xx] 
converge in distribution for A — > oo. The 
numerical results indicate that such a central 
limit theorem might indeed hold. While the 
tools used here for estimating E[Xa] do not 
give information about the distribution, it is 
reasonable that the randomness in the map 
assured by Pesin theory makes the Lyapunov 
exponents in the limit A — > c» behave like 
the Lyapunov exponents of Markovian cocycle, 
for which central limit theorems are known j5lll3|. 



9) The size of the Pesin region. 

Because the averaged Lyapunov exponent of the Standard map can not be bigger than log(A/2) -I- C(2, A), the 
conjectured bound would lead to the size of the Pesin region, the set where the pointwise Lyapunov exponent 
of the Standard map is positive, has Lebesgue measure which is bounded below by 

(log(A/2) - C(A)) _ (log(A/2) - arcsinh(l/A) ^ log(2/V3)) 

(log(A/2) + C(2,A)) (log(A/2) + C7(2,A) " ^ ^ 

For example, for A — 5.42, the Pesin region would already cover more than half of the phase space. (In 
comparison, it was measured numerically that for A = 5.0, the size of the elliptic regions is less than 2 per- 
cent of the phase space Q^j). It would follow from formula (|lUf) that the set, where the Lyapunov exponent 
is zero has measure which is smaller than 0(log(4/3)/ log(A)) for A — > oo. If C(A) could be replaced by 
c(A) = m(2. A) = log||dr\cos(a;,2/)|| dxdy = 0(1/A) the complement of the Pesin region would be the order 
0(l/(Alog(A))). 

Empirically, the set with stable behavior has been measured in |19) . Empirical formulas have to be taken with 
a grain of salt for large A because the size of the elliptic islands is usually so small that say for A — 100 already, 
a computer can hardly resolve individual islands. We know that from formula I|1U|) that the total measure must 
then be less then 3/1000. The size of the elliptic islands of a periodic orbit with period p is expected to be of 
the order , where M = sup^^j, ||(iT''(x)|| pTO]- There are examples in Hamiltonian dynamics where it is 
possible to prove that elliptic regions cover a substantial part of the phase space |108j . For other nonergodic 
Hamiltonian flows, see The Carleson problem is the question whether there are A for which Txsin is ergodic. 

10) Vlasov-Toda deformation of Twist maps. 

Other cocycles with the same bounds on the Lyapunov exponent can be obtained with isospectral Vlasov-Toda 
deformations of the random Jacobi operators L = A + V = aT+ (at)* -|- h [7J. Assume we deform L and so 
the transfer cocycle with the first Toda flow d = a{b{T) — b),b = 2a^ — 2a^(T^^), where 

a{x,y) = hi2{x,y) ^ dxdyh{x,y) 

b{x,y) = hii{x,y) + h22{T-Hx,y))^dM{x,y) + dydyh{T-\x,y)) . 

It is not clear whether the deformed operator Lt = atT + (ajr)* + bt is the Hessian of a critical point defining 
a twist map. The question is whether a generating function h can be derived from the deformed a, 6 such that 
h generates a twist map T satisfying discrete time Euler-Lagrange equations h2(T~^) + hi = 0. Because the 
deformation keeps the Lyapunov exponents of the cocycle constant, the lower bound on the entropy would hold 
also for the deformed twist maps. In any case, even if the 'isospectral deformation of the twist maps' should not 
exist in general, Toda deformation gives inflnite dimensional manifolds of cocycles over a fixed T y obtained 
from the cocycle AT,Asin of the Standard map for which one has the same Lyapunov exponent. 
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11) Continuity of the entropy. 

The map A i-^ /i(TAsin) is upper-semicontinuous. We don't know whether it is continuous. On the other hand, 
the topological entropy depends continuously on A because it is always a continuous function on C°° difFeomor- 
phisms on two-dimensional manifolds jllO| . 

M(^Ta,Acos) 

The topological entropy of surface diffeomorphisms changes 
only through homoclinic bifurcations jl51j . Motivated from 
results on unimodal maps .31) and consistent with experi- 
ments, it is reasonable to ask whether the topological en- 
tropy or the metric entropy with respect to the invariant 
Lebesgue measure depends monotonically on A. The Figure 
to the left shows the numerically computed Kolmogorov- 
Sinai entropy for 30 values of A in the interval [5.0, 10.0] 
(Lyapunov exponents averaged over 40 x 40 orbits each of 
length 10^.) 

The dependence of the entropy on A even appears to be smooth. Oscillations indeed become small for larger 
A because the established lower bound is realanalytic in A and the effective value of the entropy is in a small 
corridor log(A/2) + [— c(A), c(A)] which becomes increasingly narrow for \ ^ oo. 

12) Nonuniform hyperbolicity and homoclinic tangencies. 

Assume A is such that T\ has positive Lyapunov exponents on a set of positive Lebesgue measure. The Pesin 
region of positive Lebesgue measure is contained in the closure of transverse homoclinic points or in the closure 
of hyperbolic periodic orbits with transverse homoclinic intersections (see j65)'). We don't expect however in 
the case of the Standard map that for some A, we have uniform hyperbolicity on a set Y of positive Lebesgue 
measure. This is supported by the fact that for large A, there exists a dense set of parameter values A for which 
homoclinic bifurcations |114| happen in the Standard map family Uniform hyperbolicity on is excluded 
because the Chirikov Standard map can not be Anosov on the whole torus T^: it would be topologically con- 
jugated to a hyperbolic automorphism on the torus |96j and by homotopy to {x,y) i— > (2a; — which is not 
hyperbolic. Equivalent to uniform hyperbolicity on a set F C is the property that there exists an interval / 
containing E — Q such that / is disjoint from in the spectrum of L{x, y) for almost all {x, y) GY C T^. This is 
equivalent to the fact that the random operator L associated to an invariant ergodic set of positive Lebesgue 
measure is invertible in the corresponding crossed product algebra. It is an open question whether hyperbolic 
sets always have measure zero or one in the Chirikov Standard map case. This property is Baire generic in the 
class of measure preserving C^-diffeomorphisms 95 and probably always holds for realanalytic diffeomorphisms 
on the torus. Having excluded uniform hyperbolicity almost everywhere, the positive entropy result could pro- 
vide a new proof that homoclinic tangencies and consequently elliptic islands occur in the Chirikov Standard 
map for a dense set of parameters in [Ag, cxd). 

13) Coexistence. 

Conjecture II .21 establishes the existence of explicit real analytic maps on the torus, where coexistence of elliptic 
islands and positive metric entropy holds. Previously known were piecewise smooth maps |146l I147| or C°° 
maps |119| . Whether true coexistence (in the sense that the Pesin region as well as its complement are dense 
|137j ) can hold on open sets of the phase space for some Standard maps is not known. According to jl37j . this 
coexistence problem was posed already before 1969 by Sinai. It was measured numerically in jl41| that the 
closure of some orbits might have fractional box counting dimension. In higher dimensions, Mather conjectured 
|100 | that generically one should have transitivity in the sense that there exists orbits which are dense in the 
phase space. No mathematical results about these questions are yet available. It has also been conjectured 
that transitive components in the complement of quasiperiodic sets are trellises ii34j. closed sets obtained by 
closing the unstable manifold of some hyperbolic periodic point. For the Standard map, where we know now 
the existence of ergodic sets of positive measure, one can ask about the nature of these ergodic sets (i.e. the 
Hausdorff dimension of the closure) . Careful measurements done in |101| indicate that some ergodic components 
are obtained as the closure of some aperiodic orbits. 
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14) The complement of the Pesin region. 

An interesting question is what is left on the complement of the union of the almost periodic KAM set and the 
Pesin region. Is the dynamics weakly mixing but not mixing on a set of positive measure? While Pesin theory 
together with a genericity result for shift invariant measures |77) provide many invariant measures which have 
this property, they are in general supported on sets of zero Lebesgue measure. This observation by jSj can be 
strengthened by estimating the Hausdorff dimension of the measure. It was first observed in jl8| that there are 
invariant measures ^\ whose Hausdorff dimension goes to 2 for A ^ cx). This observation is based on Young's 
formula |153j for the Hausdorff dimension H{ii) of ji: 

H{fi) = 2h^{T)/X{fi) , 

where A(/i) is the Lyapunov exponent integrated with respect to the measure /i and where h^{T) is the metric 
entropy of T with respect to the measure /i. 

While periodic orbits are dense on the Pesin region and on the closure of KAM orbits (in the sense that for 
X in the KAM set, there exists periodic yn ^ x), the question whether periodic orbits are dense in T^ is open 
and related to the problem to analyze the dynamics on the complement of the union of the Pesin and KAM 
regions. 

15) A variational problem. 
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For fixed A, is there (a realanalytic) T e A" (rsp. y) for 
which the Lyapunov exponent of the fixed matrix valued 
map A\cos is maximal? If T is the identity, we have 
m(^Acos,t) > log|'^/2|. There are many maps RaTR-a 
with Ra(x,y) — {x + a,y) for which the Lyapunov 
exponent is bigger than log(A/2). If we fix the cocycle A\, 
how does the Lyapunov exponent depend on T^, when 
we deform r„ from the Mathieu case Tq to the Standard 
map case T\? The graph to the left shows a numerical 
computation in the case A = 3.0. The averaged Lyapunov 
exponent seems to decrease slightly when k is increasing. 



The fact that the Lyapunov exponent is close to the one in the Mathieu case k = is compatible with the 
empirical fact of being close to Aubry duality (see 3) 25)). 
It would be good to know the Taylor expansion of 

Ki-^ / \og{[A'^ .^T^{w){x,y)]ii) dxdy 

at K = 0, for \w\ < 1 which exists and converges for large A. This requires to understand the dynamics for 
complex K, where {xmUn) become complex too. 

16) Thermodynamic limit. 

For classes of symplectic coupled map lattices like finite dimensional versions of we could obtain lower 
bounds for the entropy which are independent of the dimension of the map. In the thermodynamic limit, these 
maps define homeomorphisms T on the compact metric space (T^)^ which preserve the product measure. While 
one has no more an Oseledec theorem in this infinite dimensional situation (the results in [1261193) do not apply), 
the averaged maximal Lyapunov exponent lim„^oo 'n~^ I(t^)z ^^g ||(iT"(x|| dx is defined and positive for large 
A and small e. The cocycle A = dT is now a bounded linear operator on l^(Z, C^). By Vesantini's theorem 
|121j . the map w i-^ n^^ log(p(A"(w)) is subharmonic if p{A) denotes the spectral radius of A. The estimates 
can be done directly for the operator-valued cocycle. It is not clear however whether the metric entropy of the 
homeomorphism T with invariant product measure is the limit of the metric entropies in the finite-dimensional 



28 



situations. 



17) Response formulas. 

The subharmonic estimates could apply for classes of dissipative Standard maps on the cylinder. We expect 
the averaged Lyapunov exponent to depend analytically on the map in open sets of realanalytic maps because 
our results imply that smooth observables see a pretty uniform hyperbolic behavior of the dynamics. Indeed, 
by taking periodic Rohlin approximations to the map T, we get 

/ log lldTj^.iJI dxdy > log(A/2) - C„(A) , 

where C„(A) = C(A) + (log(A/2) + C{X))/n. The hyperbolicity appears to be uniform when observing the 
dynamics with smooth observables. 

Could it be possible even that Ruelle's response formula 128 holds for some Standard maps with positive 
Lyapunov exponents? This formula says that for an infinitesimal perturbation T + ST of an Anosov map T, the 
natural T-invariant measure px will change to pr+ST = Pt + Spr with Spri^) = — X)^o PtH^ ° T") ■ div(X)), 
where X = 5T o is the vector field associated to the change of the map and where $ G C(T'') is an 
observable. The response formula would be expected to hold only for real analytic perturbations of T in A" and 
for a subset of smooth observables $ which do nowhere vanish on the torus. One has to see this question in the 
more general context of the 'chaotic hypothesis' |42j . 

18) Symmetry. 

The symmetry T{—x,—y) = —T{x,y),DT{x,y) = DT{—x,—y) which holds in the case of the Chirikov Stan- 
dard map seems to make the Lyapunov exponent extremal with respect to some parameter perturbations. This 
is supported by some numerical experiments. It is not known whether the Lyapunov exponent is continuous as 
a function of w for example when w i— s- iu, is the random operator Q , so that symmetry does not imply a local 
extremum. 
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Computing the Lyapunov exponent in two cases 
with lack of symmetry. To the right, we use 
Ta[x,y) = {2x + 3sin(3::) - y + a,x) £ X. The 
estimate known for e = a = extend. 



In the Figure to the left, we see two experiments. 
The left graph shows the Lyapunov exponent of 
the deformed cocycle 4> i-^ A{x + (p), where A{x) 
is the Jacobean cocycle of the Standard map. 
This is equivalent to a deformation of the map 
T by conjugation = R^^ o TR^. In the right 
graph, the map T is moved along a path a i— s- 
Ta = Ra o T in A", where Ra{x,y) = {x + ot,y). 
For (/) = a = 0, we have the Chirikov Stan- 
dard map case with A = 3. In the graph to 
the left, Herman subharmonicity argument shows 
that /p/x(At,3 cos (• + «)) da-log(3/2) > 0. Nev- 
ertheless, in both experiments, we never see a 
Lyapunov exponent below log(A/2). 



19) Diffusion and Sinai's (HI) conjecture. 

A discrete Legendre transform brings the Standard map into the Hamiltonian form T : (x, y) i-^ {x -\- y + 
f{x),y + f{x)) which is a map on the cylinder, the cotangent bundle T*T = T x R of the circle. Let A C 
be a measurable T-invariant set of positive Lebesgue measure. If {xj,yj) = T^(x,y) is an orbit in the 
universal cover R^'^, then Xj = yj — yj-i = f{xj) = ip{xj,yj) are random variables on the probability space 
equipped with the normalized Lebesgue measure. They have the mean E[Xj] = /rp2 fi^j) dxdy — and the 
variance Var[Xj] = /rpa ■ip{x,y)'^ dxdy. Interesting is the growth rate of = {Yl^Zo -^jY — iVn — VoY ■ Using 
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translational invariance E[XjX/] = E[Xj.4.i, the variance of Sn is 

n—1 n—1 n—1 

Var[^„] = E[(^ X,f] = nVar[X] + ^(n - j)E[XoX,] = nVar[X] + ^(n - j)Av(j) : 

where /x^ is the spectral measure of ?/' G L'^{^)i''P{x,y) — fix) with respect to the unitary Koopman operator 
3 1-^ (7(r), where r is the map r induced on fi. The Fourier transform of is/t^(j) = E[XoXj] = Cov[Xo, Xj], 
a correlation function. Let /3 be the infimum over all real numbers for which limsup„^gc n^^ 
is finite. If /? = 1, then Sn behaves like a random walk (a case, where the Xj are independent) and D — 
/x2 V'{x)'^ dx + limsup„^oc X]J=i(^ ~ i)Av'0) the diffusion constant. 

It is a conjecture of Sinai (HI) on page 144 in that there exists a set of positive Lebesgue measure for 
which /? = 1 if A is larger then Xcrit, where the last homotopically nontrivial KAM torus disappears. If fi^in) 
would decay fast enough, then D — J!p2 V'{x)'^ dx + '^JLo f^ipiJ)- First numerical experiments were done by 
Chirikov and Hizanidis |122| . Numerically, the Standard is reported to show such Brownian diffusive behavior 
for large A [861 . The fact of having positive Lyapunov exponents on a set of positive measure makes it 
plausible that the random variables Xg^Xn get decorrelated for n — » oo. When computing numerically the 
first few hundred Fourier coefficients for smaller g (like g = 5) and checking with the Wiener theorem, we got 
the impression that /iy still has some atoms (which prevents decorrelation) if = T^. Indeed, the presence 
of elliptic islands could be responsible for an almost periodic component in the Fourier transform oi fif. In 
numerical experiments, one can get rid of this discrete part of the spectrum by adding stochastic noise |85j . 
In any case, the (HI) conjecture of Sinai would be settled for A > Aq if one could show a fast enough decay of 
correlation of the spectral measure of ipix, y) = f{x) on a mixing component of the Pesin set fi. No exponential 
decay is necessary. It is enough to establish a power law decay of correlation /l^(j) = 0{j~^). This finite 
differentiability condition for the spectral measure fi^ is reasonable since the dynamics on ft is conjugated to 
a Markov chain. We even expect many spectral measures to be realanalytic leading to exponential decay of 
correlations. 

20) Dissipative Standard maps. 

The Lyapunov exponent of dissipative Standard maps numerically often satisfy the lower bound of the con- 
servative case. The Lyapunov exponent can drop however to zero. In order that the proof carries over to the 
dissipative case, it appears however, that the random variables Xj need a smooth distribution fi sufficiently close 
to the uniform distribution. Results for Henon maps lead to the expectation that for many parameters, there 
exist invariant SRB measures for dissipative Standard maps like 

Tb : {x,y) (x + y + X siii{x), b{y + A sin(a;))) , 

with 6 < 1. Results like in |15[ I103| support this. Viana conjectured that in general a map with nonzero 
Lyapunov exponents almost everywhere in the phase space has an SRB measure |142| . If this conjecture is true, 
it indicates that for most values of b, an SRB measure should exist. 

The estimates of Lyapunov exponents using spectral methods could even carry over to the case b — 0, where 
on each invariant set y — a we get the one-dimensional Arnold family a; i— > x + a + Asin(a;). If such a circle 
map has a smooth invariant measure sufficiently close to the Lebesgue measure, the Lyapunov exponent with 
respect to this measure should be > log(A/2) — C(A) for large A. 

21) More general stability of positive metric entropy? 

The stability, we have established holds for estimates of the entropy which are done using subharmonicity rsp. 
the Jensen formula. This result provokes the question, whether, in general, positive metric entropy is an open 
property in the realanalytic category: is it true that for any realanalytic, measure preserving diffeomorphism T 
on the torus with positive metric entropy, there exists a Banach space of realanalytic measure-preserving maps 
such that an open neighborhood of T has positive metric entropy? 

A question related to this stability problem is whether the Riesz measure dk of the w-parameterization of the 
operator L for a general twist map has the property that the potential / log — w'| dk{w') does not fluctuate 
too much around its mean on the unit circle {\w\ = 1}. This could be used to establish the stability if the 
entropy is comparable to upper bounds of the entropy like in the Standard map case. 
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22) The Herman spectrum. 

One can also look at different analytic paranieterizations of the cocycle. The Herman spectrum of a cocycle A is 
the set of complex numbers z, such that is not uniformly hyperbolic, where = ( z-z^'^ z+z^^ I ^ 

\ 2i 2 / 

SL{2, C) which has the property that _R(e'") — R{a) E 5*0(2, R). It is a subset of {|z| — 1} because for \z\ ^ 1, 
one can find a strict coinvariant cone field in CP^ There is a measure /i supported by the Herman spectrum 
and a harmonic function g such that the Lyapunov exponent /^(z) = /i(zi?(z)A) satisfies 

= y" log \z - z'\ d^x{z) + g{z) . 

This abstract Thouless formula [56) follows from Riesz theorem and fi is the Riesz measure of the subharmonic 
function ^{z). The proof of Conjecture ?? shows log(/i(/3)) > log(cos(/3)A/2) — 0(1/A) for all f3. The Lyapunov 
exponent is realanalytic and positive outside the Herman spectrum. At such points, the directional derivative 
with respect to variations of the angle /3 in z = re^^ can be computed using a formula of Ruelle |124j for 
the Frechet derivative of Lyapunov exponents on the open set of uniformly hyperbolic cocycles. One gets 
■jp^{f3) ~ Jrj,2 cot{uj{x,y)) dm{x), where lj{x) is the angle between the stable and unstable directions ,m~ 
at the point {x^y) (see [72]). This formula shows that the Lyapunov exponent can change a lot if the stable 

/ c b(x) \ 

and unstable manifolds are close. In the case A{x) = ( q -i ) with constant |c| ^ 1, one obtains 



d f c 

—fj,{A{f3))= cot{uj{x,y)) dxdy ^ 

d(3 J^2 1-1 



b dxdy . 



If A is the cocycle of a Standard map in Hamiltonian form like in 19), there is a spectral gap containing 

/ 2^1/2 \ 

[3 — — 7r/4 in the Herman spectrum, because R{~T:/4)dT\ ~ 0-1/2 01/2 \ t \ 01/2 )• We compute 

\ 2 ' -t- 2 ' A • cos(x) 2 ' / 

■j^A{—Ti/A) = 1. The point z = is in the Herman spectrum because it follows from [971 1127| (see also [72]) 
that z = 1 is outside the Herman spectrum if and only if the map Tx f is Anosov. The rotation number of Ruelle 
|127j defined for lifts of A into the universal cover of SL(2, R) plays the role of the integrated density of states in 
the case of the Schroedinger spectrum. It is uniquely defined up to a multiple of 27r if one fixes p{A{0)) = 0. The 
rotation number is constant on an interval / if and only if / is not in the Herman spectrum |72| . Furthermore, 
p{P) = dfi{a), showing that /i is an 'integrated density of states'. From the stable and unstable direction 
fields, one can construct (nonselfadjoint) random Jacoby operators L{z) having those direction fields m^(z) as 
Titchmarsh-Weyl functions (with energy E — 0). They are analytic in z outside the Herman spectrum. Call 
a subset of the unit circle a part of the 'absolutely continuous spectrum' of L{z) if there, m"'"(z) = m~{z). 
This absolutely continuous spectrum is a subset of the Herman spectrum. An adaption of Kotani theory |23| 
shows that the 'absolutely continuous spectrum' is the essential closure of the set, where the Lyapunov exponent 
is zero and that in the ergodic case, the existence of some absolutely continuous spectrum implies that T is 
'deterministic' (see ^njEHUHO])- We know therefore for the Standard map that there is no absolutely continuous 
Herman spectrum for almost all {x,y) in the Pesin region if A is large. We call z an 'eigenvalue' of L{z){x,y) 
if L{z){x,y) has an eigenvalue E = 0. The set of these 'eigenvalues' forms a 'discrete spectrum' of L{z){x,y) 
which is a subset of the Herman spectrum and which we expect to be nonempty for Lebesgue almost all {x, y) 
in the Pesin region. 

24) Quasiconformality. 

The map 



z^Gniz)^ / \og\z^Al{z,y)\dy 

is not conformal for small \z\ in general. For realanalytic T near a linear automorphism [x, y) i— > (a; + y, y) it is 
quasiconformal for all n. Question: is there an open set of realanalytic, measure-preserving maps T for which 
Gn are quasiconformal? If yes, it would be important to estimate the complex dilation k{z) = dGn{z) / dGn{z) 
and to estimate the Lyapunov exponent in terms of k. For k = 0, one has conformality and the subharmonic 
estimates of Herman apply. 
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25) The Aubry duality transform. 

The Aubry duality transform is an involutive map on a class of operators. It provides in the Mathieu case 
an elegant way to estimate the Lyapunov exponent. The transform can be defined in more general situations: 
Let Tfc be a cyclic interval exchange transformation on T and let -La cos, be the corresponding Schrodinger 
operator. The density of states of the random operator L is the same as the density of states of the operator 

(Lu)„(0) = u„+i(0) + u.a-i{e) + Acos(T"0)u„(0) 

on the Hilbert space H = L^(T x Z). There is a piecewise smooth potential V such that cos{T"{9)) = V{d + na), 
where a — 1/k. Let V{9) = J^n ^ exp{in6) be its Fourier series. The duality transform (see [SU)') 

satisfies tU = {X/2)Ua, where tu„(6') — Un+i{d) and aun{9) = exp(i(6' + na))u„(6'). The operator L can 
be written as L = r + r* + A K^" so that U*LU = a + a* + X/^J^n ^nir" + (t")*). We are interested 
in the density of states dk, of L, = {2/X)U*LU = e(cr + a*) + J2n'"n{T" + r*) and /(e) = log(det(L,)) = 
/log|£^| dkf{E) which satisfies /(O) ~ 0. Because log(det(LAcos,Tfc)) > log(A/2) + /(2/A), we want to 
estimate /(e) from below for small e. If Tk{x) — x + 1/fc, where /(e) is the Lyapunov exponent of a 
symplectic transfer cocycle, we have /log|i?| dk^{E) > 0. In general, the Thouless formula just becomes 

log(det(Le)) = limjv_»oo /t ^°s('^j (■^' ^)) '^^^ where Aj(e,a;) are the eigenvalues of the truncated N x N 

matrix Le^Af (a;). While /(O) = and a perturbation lemma of Lidskii ( [1291 1 assures that |Aj(e,a;) — Aj(e,0)| < e, 
this does not allow us to estimate /(e) from below. The perturbation problem to estimate /(e) seems still diffi- 
cult and it is not clear, whether the Aubry transform, which transformed the perturbation problem V i— s- V + eA 
into a perturbation A i— > A + eV has made things simpler. 

26) Lower bound for topological entropy? 

The topological entropy of T g A" on is bounded below by log(sp(r*)), where sp(T*) is the spectral radius 
of : i7*(T^) I— > _ff*(T^) |1U2) . Can the metric entropy of a (smooth) T £ X with respect to the invariant 
Lebesgue measure become smaller than sp(r*)? 

27) Conjectures on the Standard map. 

The following of conjectures about Standard maps have still to be settled. 

I) In Pni, the entropy of the Standard map Tasiii was measured > log(A/2). Chirikov formulated the possibility 
that the elliptic islands might cover arbitrary large regions of the phase space. 



I J72|: The Kolmogorov-Sinai entropy of the Chirikov Standard map is 
> log(A/2). 



II) |134| introduced the problem of determining the spectrum of random operators {L{x, y)u)n = Un+i — 2w„ + 
Un-i + A cos(xn)w„, where T"{x,y) = {xn{x,y),yn{x,y)) is the orbit starting at {x,y). 



II [1 34| : The operator L{x,y) has some point spectrum for a set {x,y) of 
positive Lebesgue measure if A is large enough. 



Ill) '16' asked whether there exists a set of parameters A with full density at oo for which the Chirikov Stan- 
dard map has no elliptic islands. We know that we have for large A an open dense set of parameters with no 
ergodicity and positive entropy. 
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Ill [16 : There are parameters A with full density at oo for which the Chirikov 
Standard map Txsin is ergodic. 



Heuristic arguments in 07] predict that the Lebesgue measure of the set of parameters AG [r, r + 1] which lead 
to nonergodic Standard maps TAsin is of the order 0(l/r). 

IV) The following problem is related to III): 



IV | 132 |: For large A, the Standard map has (in some Baire topology of 
maps) a neighborhood in which a residual set of /'s gives ergodic maps T/ . 



V) In the textbook |131| p. 144, conjecture (112) contained as a second part the statement that the entropy grows 
to infinity for A — > oo. While we solved this part of the problem here, the first part of the conjecture is still open: 



V [131]: the entropy of the Chirikov Standard map is positive for all A > 0. 



One should be able to prove that for any real-analytic non-constant periodic /, there exists Ao > such that 
the Standard map Txf has positive Kolmogorov-Sinai entropy for all A > Ao and that for any real-analytic, non- 
constant, periodic map / : T'' and every symmetric, constant matrix E e GL{N,Z), all the Lyapunov 
exponents of the symplectic map TEx+\f are nonzero on a set of positive Lebesgue measure for large enough 
|A|. All these properties are expected to be stable with respect to realanalytic perturbations of the map Tex+x/, 
and hold for a fixed cocycle for all T G 3^. 
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